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In a recent work Povolotskj® provided a three-parameter family of stochastic particle systems 
with zero-range interactions in one dimension which are integrable by coordinate Bethe ansatz. 

Using these results we obtain the corresponding condition for integrability of a class of directed 
polymer models with random weights on the square lattice. Analyzing the solutions we find, besides 
known cases, a new two-parameter family of integrable DP model, which we call the Inverse-Beta 
polymer , and provide its Bethe ansatz solution. 

PACS numbers: 


I. INTRODUCTION AND MAIN RESULTS 
A. overview 

There is considerable recent interest in exact solutions for models in the universality class of the ID stochastic 
growth Kardar-Parisi-Zhang equation (KPZ) 1 . Models in the KPZ class share the same large time statistics, also 
found to be related to the universal statistics of large random matrice:P. Methods developped in the context of 
quantum integrability are exploited and broadly extended to solve a vari ety of ID stochastic models. The Bethe 
ansatz solution of the attractive delta Bose gas (the Lieb-Liniger model®!) was combined with the replica methocP, 
to obtain exact solutions for the KPZ equation directly in the continuum and at arbitrary time, for the main classes 
of initial conditions (droplet, flat, stationary, half-space(P®l The Cole-Hopf mapping h ~ In Z is used, where h 
is the height of the KPZ interface and Z the partition sum of a directed polymer in a random potential (DP). 
Hence in the continuum, studying KPZ growth is equivalent to studying the DP model, an equilibrium statistical 
mechanics problem with quenched disorder. The time in KPZ growth becomes the length of the polymer t. The 
replica Bethe ansatz (RBA) method then allows to calculate the integer moments Z n and, from them, to retrieve the 
probability distribution function (PDF) of Z. Since the last step is non-rigorous because of the fast growth of these 
moments, the mathematical community has concentrated on the exact solution of discrete models, which in favorable 
cases, do not suffer from the moment growth problem. Discrete models, such as the PNG growth modeP^HUl, the 
TASEP and ASEP particle transport riiodeP®Sl anc j discrete DP modeld 17 21 ! 22 ! played a pionneering role in unveiling 
the universal statistics of the KPZ class at large time (the Airy processes). Recently they have been considerably 
generalized, unveiling a very rich underlying ’’stochastic integrability” structure 23 33 ! Since in suitable limits (e.g. 
ASEP with weak asymmetry, q -TASEP with q —i 1, semi-discrete DP) they converge to the continuum KPZ equation, 
they also led to some recent rigorous results for KPZ at arbitrary tim<P®®! 

Besides their interest in relation to KPZ growth, directed polymers are also important in a variety of fields. This 
includes optimization and glasseiP® 1 !, vortex lines in superconductor^^, domain walls in magnetJ®! disordered 
conductors®!, Burgers equation in fluid niechanicPP exploration-exploitation tradeoff in population dynamics and 
economics— and in biophysics 4 ®®. In some situations (heavy tailed disorder) they exhibit anomalous (non-standard 
KPZ) scaling^® 51 . Apart from models on trees, exactly solvable models of DP (e.g. on regular lattices) remain, 
however, exceedingly rare. We will present in this paper a new solvable DP model. 

On the square lattice a few remarkable solvable DP models have been found. The first that was discovered is at 
zero temperature T = 0 (i.e. it amounts to find the minimal energy path, energies being additive along a path), with 
a geometric distribution (of parameter q < 1) of on-site random potential^®. The second that was discovered, called 
the log-Gamma polymer 5 ?!, j s a finite temperature model (as it focuses on Boltzman weights, which are multiplicative 
along a path), with a so-called inverse gamma distribution for the on-site random weights, with parameter 7 . This 
weight distribution has the peculiarity of exhibiting a fat tail P(w ) ~ w~ 1+1 . These models are not unrelated: in 
the limit 7 —> 0 (so-called zero temperature) the lo g-G amma converges to the q 1 limit of the Johansson model 
(i.e. with exponentially distributed on-site weights)®!. They were both proved to belong to the KPZ class, with 
convergence of the free energy PDF to the GUE Tracy Widom distribution. The Johansson model was solved as a 
determinantal procesd®!. The log-Gamma model was solved using the gRSK correspondence (a geometric lifting of 
the Robinson-Schensted-Knuth (RSK) correspondence) leading finally to an expression for the Laplace transform of 
P(Z) as a Fredholm determinant ! 53 ! 54 !. 

Recently, we provided a solution of the log-Gamma polymer using replicas and the coordinate Bethe ansatz, closer 
in spirit to the integrability methods used to solve KPZ 55 . As in the continuum, this replica Bethe ansatz approach 
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consists in computing the moments of the partition sum of the DP using a transfer matrix (i.e. recursive) formulation of 
the problem. This formulation can formally be interpreted as a discrete-time quantum mechanical model of interacting 
Bosons. Such a connection between discrete-time particle models and lattice DP was also noted, and exploited in® to 
unveil and study a new integrable DP model with Gamma distributed Boltzmann weights, called the “Strict-Weak” 
DP model, as the q —> 1 limit of the discrete time q-TASEP modef®. In parallel, this DP model was also solved using 
the gRSK correspondence irf®. 

In a recent seminal work, Povolotskj® provided a three parameter family of discrete-time stochastic interacting 
particles systems with zero range interactions (“zero-range processes” (ZRP)) called the (q, /z, i/)-Boson process and 
integrable by coordinate Bethe ansatz. This led to further rigorous work on this class of particle model and on a dual 
model, termed the q -Hahn T ASEP, which eventually allowed to unify integrability properties of ASEP and q-TASEP, 
a long-standing goa l 25 * 30 31 . On the directed polymer side, this work also led to the discovery of a new integrable 
model, called the “Beta” polymer, introduced and studied in®. There the model was solved as a q —> 1 limit of the 
(q, fi, i/)-Boson (in analogy with the Strict-Weak case), but the authors also already provided a direct replicas Bethe 
ansatz solution of the model. 

The aim of the present paper is to explore more systematically the consequences of Povolotsky’s work to directly 
search for, and attempt to classify, the corresponding family of integrable DP models. Integrability then leads to a 
constraint on the integer moments of the Boltzmann weights distributions, and we search for solutions in terms of 
PDF of bond and site disorder. 

We find that there are two main solutions, the first one corresponding to the Beta polymei^®. The second however 
is new and corresponds to weights v on horizontal bonds, and u on vertical bonds of the square lattice, with the 
following PDF: u, is distributed according to: 


PyA u ) 


r(7 + l) 1 (, _ 1 
T( 7 )r(/3) v.i+'r \ u 


u e [l,+oo[ , 7,/3 > 0 


(1) 


The weights are correlated on bonds which share a top/right site (see Fig. [I]), with v = u — 1 € [0, +oo[ but otherwise 
uncorrelated. Given the form of |l]) we call our new model the Inverse-Beta polymer 63 . 

We will provide in this paper the coordinate Bethe ansatz solution to this model, as well as some explicit integral 
representation and Fredholm determinant formulas for its Laplace transform. It is interesting to note that for /3 —> +oo 
this model, under suitable rescaling, converges to the log-Gamma polymer (see below). Hence it can be considered as 
a generalization of the log-Gamma polymer. 


B. Main results and outline of the paper 


The first result of this paper, obtained in Section [TT] are some general conditions for a finite temperature model 
of directed polymer on the square lattice to be integrable using the coordinate Bethe ansatz. The only hypothesis 
are that Boltzmann weights on horizontal edges and vertical edges can be correlated only if they share the same 
top or right site (an example of short-range correlations), and that they are homogeneously distributed. Within this 
framework, in Section m we attempt a classification of integrable DP models, retrieve the known integrable models 
and introduce a new one, the Inverse-Beta polymer, whose Boltzmann weights are distributed as ([I]). This model has 
two parameters 7 , /3 > 0 and contains the log-Gamma and Strict-Weak polymers as scaling limits. More precisely, we 
show that the partition sum Z t (x) of the Inverse-Beta model (see Section E3 for the definition) converges in law to 
the partition sum of the log-Gamma (resp. Strict-Weak) polymer Z^ G (x) (resp. Zf w (x)) as 


A" 


7 lg 


0 ) 


lim 7 : 

7—>-oo 


Or) 


( 2 ) 


In Section [TV] we use the coordinate Bethe ansatz to study the Inverse-Beta polym er with point-to-point boundary 
conditions. We obtain an exact result for the integer moments of the partition sum Z t (x) n (60), defined for n < 7 . 
Using this result, we conjecture the formula (69) that expresses the Laplace transform of Z t {x) as a Fredholm 


63 Note that a nomenclature based on the names of the weight distributions, the log-Gamma polymer could be called the Inverse-Gamma 
polymer, and the Strict-Weak the Gamma polymer. Alternatively our model could be called the log-Beta polymer. 
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determinant e~ uZt< ^> = Det (I + Ktx) with 

r+ °° dk -1 f ds 
-oo 77 2i Jc sin(Trs) 


K t , x {v i,v 2 ) = 


u s e -2ik(v 1 -v 2 )-s{vi+V2) 


(3) 


r (-f + %-ik) 


r (i 


— ik) 


1+X 


r (-| + ?+ **) 


r (i 


ik) 


1 -x-\-t 


T(/3 + ik+^ + f) 
T(/3 + ik+l~ f) 


where C = a + ® with 0 < a < min(l, 7 ) and K t x : L 2 (M + ) —> L 2 (M + ). Alternatively, we obtain an equivalent 
Fredholm determinant for the same quantity with a different kernel (which contains notably one less integral) in (|70| ). 
By analogy with a known formula for the log-Gamma polymer, we also conjecture a n-fold integral formula ( |71[ ) for 
the Laplace transform 


e —uZ t (x) — _ 



TT +7. TT 1 

(my ~=i 2in j}k=i T ( w i “ 


J r[a - WjY 


j 7 r (l + a-Wj) 

r(7) 


1 -a + P) y +J 2 j 


(4) 


with 0 < a < min(l, 7 ), valid for Re(u) > 0, 1 < J < I and where x = I — 1 and t = I + J — 2. Using an 
asymptotic analysis of our Fredholm determinant formulas, we show in Section |IV C the KPZ universality of the 
model for polymers of large length t —> oo with an arbitrary angle ip s] — 1/2,1/2[ with respect to the diagonal. More 
precisely, we show 


lim Prob 

t—> OO 


10 gZ t ((l /2 + lf)t) + tCy 

Xrn 


<23 z) = F 2 (z) 


(5) 


where F 2 (z) is the standard GUE Tracy-Widom cumulative distribution function, X v ~ + 3 and the ((^-dependent) 
constants are determined by a system of equations (801 that involves the digamma function ip. As a particular case 


we study these characteristic constants for long polymers with the ‘optimal angle’ ip = ip* (in the sense that the mean 
free energy c v is minimal for this angle) and find explicit expressions as 


ip* = 


1 ^(/g + 7/2) 

'2 i//( 7 /2) 

= +7/2) ^ i’iP d 
K* = + 7 / 2 ) - V ,// (7/2)) 


< 0 

•7/2) 


1/3 


( 6 ) 


Finally, in Section (IVD) we study a two parameters zero temperature DP model that we obtain as the limit 
7 = £ 7 ' and /3 = e/3' with e —> 0 of the Inverse-Beta polymer. This study is close in spirit to the one made iiJSH 
where the zero temperature limit of the Beta polymer is studied, but the models are qualitatively very different. 
The energy of this model are distributed as (£(,,£') = (—£)£+ (1 — ()Ey — QEy) where (£' u ,£' v ) are the energies 
on vertical and horizontal edges, £ is a Bernoulli random variable of parameter p = Q' j(p/ + j3') and Ey and Ey 
are exponential random variables of parameter 7 ' > 0 and (3' > 0, independent of </. This model generalizes the 
known zero temperature limit of the log-Gamma directed polymer and we obtain exact results for the cumulative 
distribution of the optimal energy, noted £(t, x ), of this zero temperature model Prob(<£^ t ,x) > r )■ I n particular we 
obtain a Fredholm determinant formula (99) Prob( <£(t,x) > r) = Det (/ + with 


r+oo 

K Zx°(v i,u 2 )= 77 


dk 
7 r 


ds 


0 sr—2ik(y\—V2)—s(v\+V2) 


I + l~ik 


c 2*7rs 

1 + X 
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where (7 = 0 + ® with 0 < a < 7 ' and Kf~° : L 2 (M + ) —>• L 2 (M_|_). We also conjecture an equivalent n-fold integral 


formula (1011 
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with 0 < a < 7 '. Using our exact results, we conclude this section by showing the KPZ universality of the zero 
temperature model in (106 1 . In the case /3 '/Y —> oo the model maps onto the Johansson DP model with an exponential 


distribution and we show that our solution reproduces all the (non trivial) angle dependent normalizing constants in 
the statement of convergence to the GUE Tracy Widom distribution. 

A series of appendices also contains additional discussions and some technical details separated from the main text 
for clarity. 


II. DIRECTED POLYMERS ON THE SQUARE LATTICE: REPLICA METHOD AND 

INTEGRABILITY 


A. Definition of the model 

We consider the square lattice Z 2 with coordinates (f, x) with x the usual horizontal coordinate and t a coordinate 64 
running through the diagonal of Z 2 as depicted in Fig. [l] We will also sometime use the usual euclidean coordinates 
(/, J) on Z 2 with x = I — 1 and t = I + J — 2. The first quadrant is thus /, J > 1 and x, t > 0. A directed polymer 
model on Z 2 is defined by the partition sum 

Zt{x)= n-e (9) 

7r:(0,0)—qt,x) ee7T 

where the sum is on all directed (i.e. up/right) paths with fixed starting point (0, 0) (corresponding to (/, J) = (1,1)) 
and endpoint ( t,x ), and the product is on all edges e = (t\ x') -A (t' + l,x') or e = (t',x r ) -A (t’ + \,x’ + 1) visited by 7r. 
Here for definiteness we consider a directed polymer model with fixed endpoints, but the model can be generalized to 
other boundary conditions. We also restrict ourselves to models with on-links Boltzmann weights w e . Obviously, by 
redefining the weights w e , one can also include on-sites Boltzmann weights so that this hypothesis is non restrictive. 
The Boltzmann weights are positive random variables w e € R + . We will generally note u (resp. v) the Boltzmann 
weights on vertical (resp. horizontal) edges: 

w e = u t , x if e = (t — 1, a;) -A (f, x ) , 

w e = v t , x if e = (t — 1, x — 1) -A (t, x) . (10) 

We will consider the class of models with the following structure of local correlations, which naturally emerges in the 
integrable family we are studying: the weights on edges arriving at different sites are statistically independent, but the 
weights of two edges arriving at the same site are correlated. Thus the model is defined by a (common) joint PDF for 
the weights of the type ( Ut iX ,Vt iX ), denoted p(u,v) (see Fig. [lj. The pairs (ut tX >Vt, x ) are chosen independently from 
site to site. In the following, the overline (.) denotes the average of a quantity over all realizations for the ( Ut, x , v t,x)- 


B. The replica method and the coordinate Bethe Ansatz. 

In general, one is interested in computing the PDF of Z t (x) or of its logarithm \ogZ t {x). The replica method 
consists in first studying the equal-time moments of Z t (x): for n £ N and whenever they exist, one defines 


• • • ,%n) = Z t (x i)... Z t (x n ) . (11) 

In the general case these moments are only defined for n < n max , because of possible fat-tail in the distribution of 
Boltzmann weights, such that u ni v n2 are finite for ni+ri 2 < n max but are infinite for some (ni, n 2 ) with ni +712 > n max . 
In the log-Gamma case one has n max = (7 _ J whereas in the Beta polymer case n max = 00 . Obtaining the PDF of 
Z t (x) from the knowledge of the moments is usually non-trivial, especially when n max < 00 . In this case the procedure 
is non-rigorous and one has to perform some analytical continuation as in the log-Gamma case. In most of this section 


64 Note that the (t.,x) coordinates of the present paper do note coincide with the ones offSthat we denote (T,X). To compare formulas, 
one can use t = T and x = t/2 + X. 







5 



FIG. 1: General scheme for the models of directed polymer studied here. Blue (resp. Red) : couple of correlated Boltzmann 
weight on edges arriving at (t = 6,x = 1) (resp. (t = 6, x = 4)). Green: two admissible (i.e. up/right) paths for polymers with 
starting point (0,0) and endpoint (8,4). 


we will not discuss this issue and only focus on computing the moments for n < ro max , which is a well-defined problem. 
The problem of computing ip t is manageable thanks to a the recursive formulation of <©: 

o(^) - ^£E,0 

^t -\-1 (*^) — 'U , t-\-l,xZt ( 3 ?) + (*£ 1 ) • ( 12 ) 

This can be translated to a recursive (i.e. transfer matrix) equation for ij) t : 


o(^l • • • 1 %n) — $xi,0 • • • ^x n ,0 

4>t+i(xi ...,x„)= Y aS xY,-^x„^t{x 1 - <5i, • • ■ , X n - 6 n ) = {T n lf>t)(x 1 ■ • •, x n ) 

{<5l,— ,<5n}e{0,l} n 

a t \= ]^[ . (13) 


Where we used the statistical independence of the Boltzmann weights ending at different sites and the definition of 
( u,v ). Note that the evolution equation (13) is symmetric by exchange xi o Xj. Therefore, since the initial condition 
is also fully symmetric, if one is able to find all the symmetric eigenfunctions ip M of T ni i.e. a complete basis of 
symmetric functions such that T n ip M the problem is essentially solved. In the already known models, it was 

possible to find the eigenfunctions of T n in the form of the coordinate Bethe Ansatz. More precisely, in the sector 
W n = {xi < • • • < x n } (this defines the Weyl chamber), one looks for eigenfunctions of the form 


i>p(xi, ...,x n ) = ip^x 1 ,... ,x n ) if X! < ••• < x n 

n 

i> ll (x 1 ,...,x n ) = y A °n z *\i) ( i4 ) 

<Tes n i=1 


where the sum runs over all permutations of {1,..., n}, and the variables A a and Zi are complex numbers. The wave 
function ip^xi ,..., x n ) is deduced in the other sectors by using that it is fully symmetric function of its arguments. As 
one can guess, this form of eigenfunction is restrictive and it can only works if the variables a ^], or equivalently 
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the integer moments u ni v n2 for ( 711 , 712 ) £ N 2 obey a particular structure. Thus one can hope to classify the models 
that are solvable by the coordinate Bethe Ansatz. In fact, (13) is reminiscent of the equations usually considered in 
the study of zero-range stochastic particle systems, for which a classification was proposed in^S and latter extended 
irP^. In the next section we follow the route oP^ and adapt it to our setting to deduce a classification of integrable 
directed polymers models 65 . 


C. The constraint of integrability on integer moments. 


If one can diagonalize the evolution equation (13) using the Bethe Ansatz ( |14| ), then in the sector Wn= {%i < ■ ■ ■ < 

1 particles sit on distinct sites and do not interact), one must have 


x n } (i.e. the interior of the Weyl chamber where a 

)X\ j ■ ■ • j X n ) — (Tri^Pfj) (ccr , . . . , X n ) 

{<5i,...,<5 n }e{0,t}'‘ 

n 

= ]J(u + 1, • ■ ■ , Xn) , 

i= 1 


(15) 


and this already imposes the eigenvalue to be A M = Jj" =1 (u + vz i ). Note that this is a direct consequence of the 

O 

weights having zero-range interaction: in the Wn sector, the operator T n just acts as a (biased) diffusion operator 
on the one-dimensional line. Let us now look at what happens when exactly two particles are at the same position: 
X\ < ■ ■ ■ < Xi = Xi + 1 < ■ ■ ■ < x n . In this case, the evolution equation reads 


A^t/^aq,... ,x n ) = 


^2 («) n 2 (u 2 ip^(x 1 -6 1 ,...,x i ,x i ,...,x n -S n )) (16) 

{« lt — ,<5,d—{<5i,a i+ i}e{o,i}"- 2 

+2mi/j^,(xi - Si,... ,Xi - l,Xi,... ,x n - 6 n ) + v 2 ip fJ ,(xi - Si,... ,Xi - l,Xi - 1,... ,x n - S r 

Where we used the symmetry of to express each terms with coordinates in the Weyl chamber W n . However, the 
left-hand side of (16) is already constrained to be equal to the right-hand side (last line) of (fl5| even for Xi = Xi + i 


because the eigenvalue A M is entirely determined by (15). For this equality to hold Vaq < • • • < aq = aq+i < • • • < x r 
for an eigenfunction of the form (14) one must have, rewriting (16) in terms of 7 /^, 


« 2 Vv(z 1 - 61,..., Xi, Xi,..., x n - S n ) + 2 uv'ip l j i (xi - ... ,Xi - l,x h ... ,x n - 5 n ) 

+v 2 p) fJ ,{x 1 - <5i,..., Xi - 1, Xi - 1,..., x n - 5 n ) 

= {v) 2 p)^{xi - <Ji,.. .,Xi - l,Xi - 1,... ,x n - 5 n ) + (u)(w)V>ft(®i - 5i,... ,Xi - 1 ,Xi,.. 

—\ 2 „ 


-s n ) 


+{u)(v)ip^x 1 - 61 ,... ,Xi,Xi - 1,... ,x n - 8 n ) + (uyip^X! - Si,... ,Xi,Xi,... ,x n - 5 n ) . 


(17) 


Notice that the third term in the right-hand side in involves coordinates outside of the Weyl chamber and is thus 
not a physical term. For the two particles problem to be solved, it must have the value 


Vv(£j, Xi - 1) = a ip^Xi, x^ + b i’^Xi - 1, x^ + cip^Xi - 1, aq - 1) 

■-(u) 2 

—\ /—\ 5 


a = 


u 2 — (u ) 2 2 uv — (u)(v) 

b = 


{u){v) 


(u)(v) 


c = 


(u)(v) 


(18) 


where here, for clarity, we did not write the other particles positions. For obvious reasons, (18) is called the two- 


particles boundary condition. It is also a consequence of the short-ranged nature of the correlations between Boltzmann 
weights that the two-particles evolution equation can simply be interpreted as a two particles boundary conditions. In 


65 Note that the extension ofMJ Q f t.h e classification o^l corresponds to stochastic particles systems with non-simultaneous updates. Hence 
the classification off^is sufficient in our context. 
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terms of the Bethe Ansatz (14), it imposes the quotient of two amplitudes A a related to each other by a permutation 
to be given by 


S{Zi,Zj) := 


A 




c + bzj + aziZj — Zj, 


A...ij... c + b Zi + aziZj — Zj 
Where this defines the S matrix. This can be solved as 

C T ^~tx(*) ^Ax(i)Ax(j) Za(j) 


Ax = e(a) 

l<z<j<n 


c + b Zi + aziZj — Zj 


(19) 


( 20 ) 


As a consequence, up to a multiplicative factor, the form of the Bethe Ansatz is now entirely specified and something 
special has to happen if it also solves the m particles problem (case where m particles are at the same position) 
for arbitrary 2 < m < n. Indeed, for arbitrary m, one can repeat the same analysis and check that the ansatz 
(14) works, i.e. that the evolution equation with m particles at the same position (generalization of (16)) can be 
transformed into the free evolution equation (generalization of @) by only applying ( p~8] ) recursively. Schematically, 
this is conveniently encoded in a non-commutative algebras with two generators (A, B ) such that 


BA = a A 2 + bAB + cB 2 


( 21 ) 


which encodes what happen when one transforms a forbidden term of the form .. ,Xj, Xj — 1 .. .) into a sum of 
terms with coordinates in the Weyl chamber. In this language, the model is indeed integrable with the coordinate 
Bethe Ansatz (14) if and only if 


(uA + vB) m = E u ” 1 v m ~ ni C ™ 1 A ni B n 


( 22 ) 


ni —0 


Where the right-hand side represents the true evolution equation that only contains terms in the Weyl chamber, i.e. in 
this language, that only contains ordered words of the form A ni B n2 , and the left hand side is the formal free evolution 
eq uati on, which contains various terms outside the Weyl chamber, i.e. wrongly ordered words. The right-hand side 
of (22) can be computed using the formula appearing in 24 l In the context of this paper, only models satisfying the 


“stochasticity hypothesis”, a + b + c = 1 (corresponding to a conservation of probability) were considered. Here in 
general, this hypothesis has to be relaxed and a + b + c^ 1. This is easily done by a scale transformation 66 , i.e. we 
introduce new parameters (p, a', b', c') and generators ( A',B '), such that 


a = a p 
A' 

A = — 
P 


b = b' 

B = B' 


c = 


B'A' = a'A ' 2 + b’A'B’ + c’B’ 2 


where p is chosen such that a' + b' + c' = 1, in which case we can use the results of 
generators ( A',B') and parameters (a', b', c') (called there (A,H) and a, (3, 7 ). We obtain 


(23) 

(24) 

24 in terms of the new 


(uA + vB) m =(-A' + vB') m 
P 

= (p +0 ) ^ A ' + ^ 


= -+S E ^Ani\m)(ArHBT- ni 


ni —0 

_ \ m m 

-+v ) E P ni ^Ani\m)A^B^ 

P ‘ n,=0 


(25) 


3 We thank A.M. Povolotsky for this remark. 













where we used the same notations introduced in RefP^ for the three parameters of the model g, p, v and the auxiliary 
parameter p , so that 


P = 


- q) 


u + pv 

<t>q,vA n ll m ) = M ni — 


b' = 


1 — qv 1 — qv 

(«; q)ni {P'i q)m— m (g; g) r 


9 -^ , 1 - 9 - /1 \ 

c = -- , /i = p+i/(l-p) 

1 — gt/ 


(26) 


(g;g) ni (g; q) m -n 1 

from equations (8) and (26-27) in RefpS The g-Pochhammer symbol (used extensively below) is defined as, for n > 0: 


(a;q)n = n^ 1 “ a ^) ’ (a;<?)-n = “ aq ^ 1 


(27) 


k=0 


fc=l 


and (a;g)o = 1. For given parameters q, p, v the first equation in (26) fixes the scale factor p as a function of u/v. 
Comparing this equation with (22) one sees that one must have, V(ni,7i2) such that ni + n 2 < ?i max , 


u ni v n2 = (-) ni (- 

P 1 ~P 


\ n 2 ^ / J \ n 1 ^ ,q ' )ni ^' ,q ^ n2 ( q > q ) 


n 1 +n 2 


1 


(ffi q)ni +n 2 (q;q)m{q;q) n 2 ^n±-\-n 2 


(28) 


We can now explicitly check that for ( 711 , 712 ) = (1,0), (0,1) the r.h.s gives u and v and that for ( 711 , 712 ) = 
(2, 0), (1,1), (0,2) it yields second moments compatible with all relations (18) and (23), (26). Thanks to the above 
construction based on Ref. 24 we know that it solves the integrability constraint for all higher positive integer moments 
with 71i + 712 < ft max- 

In this expression, the power-law parts are unimportant (they can be absorbed into a rescaling of the Boltzmann 
weights which cannot break the integrability of the model). We can now reverse the construction and study a polymer 
model defined with weights with moments given by 


it™ 1 v n2 = 


{q)n 1 (p'i q)n 2 (g; g) 


n 1 +n 2 


1 


( v \q)ni 


+n 2 


{q\q)nM,q)n 2 c 1 " 1 


:= 4’ q ^,u{ni,n 2 ) , 


(29) 


n 1 +n 2 


where (g, p, v) £ M 3 to obtain real Boltzmann weights. This model is automatically integrable, and with the hypothesis 
that we made, it is the only form for the moments that leads to integrability. However, we now need to check if this 
DP model really exists, namely that (29) corresponds to the moments of a PDF p(u,v). 


Let us now define the moment proble m t hat we must now solve. We are interested in finding a joint PDF p(u, v ) 
with positive integer moments given by ( |29| ) and random variables (it, v) living in one of the four quadrants (K 41 , K 14 ). 
Indeed, if that is the case we automatically find, using a change of the type ( u,v ) —> (±u, ±u), positive random 
variables with moments given by (29) (eventually multiplied by additional power laws (±l) ni (±1)™ 2 which do not 
spoil integrability). Since we extend our search to also include polymer models with 7 i max < 00, we will generally look 
for PDF with moments given by (29) for n\ + n 2 < 7i max for some 7i max - Note that if n max < 00, this replica Bethe 


ansatz method allows us to compute a-priori only a few integer moments of Z t {;x). The ultimate goal of computing 
the PDF from this knowledge is not a mathematically well-posed problem. Fortunately, as e.g. in the case of the 
log-Gamma polymer (se^ 5 for more details on this issue), and in the case of the Inverse-Beta polymer studied below 
in this paper, the situation turns out to be more favorable. Indeed in these cases, though n max < 00, the complex 
moments u ni v n2 of the PDF p(u,v) exist for (7i!,7i 2 ) in a large domain of the complex plane plane C. These are 
given by an analytical continuation of (29) and allow, using a Mellin-Barnes type contour integral formula, to recover 


the Laplace transform of p(u, v) in a rigorous manner. In this paper and as in the log-Gamma case, we adapt this 
observation to conjecture a formula for the LT of Z t (x) by using an analytical continuation of the formula for the 
integer moments that we compute using the replicas Bethe ansatz. 

The search of such a PDF p(u,v) with moments given by (29) is in general a difficult task. Notice however that it 


is sufficient to examine the case |g| < 1. Indeed, using that (i; ^)„ = (—a) n q ' 2 ( 0 ; q) n one easily sees on (29) 


that the simultaneous change q —> 1/g, p —> 1/p and v —> 1/v in ip q>ll>v {n 1 , 112 ) just multiplies it by power-law terms, 
easily absorbed in rescaling of the variables (it, v) and which cannot break the integrability of the model. 


III. INTEGRABLE POLYMER MODELS 

A. The |g| < 1 case. 


Without loss of generality we restrict ourselves in the following to |g| < 1, and we further restrict to g, p,v £ R. 
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We now consider the case where the moments exist at least up to the second moments (i.e. n max > 2). Let us 
consider the random variable z x = u + xv , x £ R. A simple calculation from (29) gives that its variance is: 


72 —2 _ (M- 1 )( 1 -4)(m-i')(/^-1)(m:e-z') 


zt - z x = 


/i 2 (i/ — 1) 2 (z^(7 — 1) 


(30) 


Under our assumptions this expression must be positive. Since the polynomial in x changes sign at x = ^ and x = ^ 
this clearly rules out the generic case for q < 1, /j,, v. 

We must thus look for degenerations with v/\x = 1 /^l so that the variance of z x can eventually be positive Vx. The 
various cases are studied systematically in Appendix [A] where we show that the only possibility for the existence of 
such a PDF is in the q —> 1 limit which we now study in details. Moreover, we also show there that the limit q —> 1 
and /x, v —> 1 at the same speed than q , contains all the interesting cases. 

Finally, note that the above considerations do not rule out completely the existence of an integrable polymer model 
with q < 1 since it could correspond to a model with n max < 2. From the discussion of the previous Section, this 


would involve however an exhaustive study of the possible analytical continuations of (29) which goes beyond the 
present work. 


B. The q —> 1 limit 


1. Form of the moments. 


Let us now discuss the q —> 1 limit. We use that at fixed n,a, ( q a \q) n —q^i (1 — <l) n {. a )m where (a) n = a(a + 
l)..(a + n — 1) is the standard Pochhammer symbol. This is easily seen setting q = e -e and taking e —>■ 0. The ratio 
( ' q ' q ' n - thus tends to the standard binomial coefficient C™. 


(g;g)m(g;q)r>-« 

To obtain a meaningful limit we scale v = q a+ P , 


F = q 0 - 


u n lv n 2 = (ei) ni (e 2 ) T 


In this case, one gets as q 

2 (oQni 

{cx T /3) ni ^_ n2 


1: 


(31) 


where we have added two power law terms with (ei, e 2 ) £ {—1, l} 2 . We were allowed to do it if we start to examine 
the moment problem for real variables. These two additional parameters are then tuned so that (rt, v ) are positive 
random variables. Since they are a-priori arbitrary we must examine all cases. Other interesting limits can also be 
considered but they can all be obtained from (31) as a new limit (see Appendix [A]) . Note that (31) implies, Vn € N, 


(u/e 1 +v/e 2 ) n = > Cl 1 


E 

ni=0 


(®)ni (^)n—r 

(a + /3)„ 


= 1 


(32) 


so that, except maybe in some marginal cases discussed in Appendix [B] this implies that u and v are correlated as 
e 2 it + e±v = eie 2 . In Appendix [Bj we initiate a more systematic study of all possible cases as are varied. 


2. The Beta Polymer and the Strict-Weak limit 


The case of (ei, e 2 ) = (1,1), a > 0 and /? > 0 indeed corresponds to the moments of two positive random variables. 
In this case, one has v = 1 — u and u £ [0,1] is distributed according to Beta random variable: 

Beta(a,P) p a}f 3 (u) = - u ) 31 

v = l -u , (a)nMn 2 _ ( 33 ) 

(a + /?) ni+n2 


Where from now on ~ means distributed as or the equivalence in probability. Note that Beta distributions satisfy 
Beta(a, f3) ~ 1 — Beta(/3, a), and that in the Beta polymer model, interverting horizontal and vertical edges amounts 
to permute a and 

Note that for this distribution of (it, v), the formula for the moments (33) can be extended to the complex moments 
and admits a more general expression as 


r(a + /3) r(a + si)r(/3 + s 2 ) 
r(a)r(/3) r(a + P + Si + s 2 ) 


( 34 ) 
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which is valid for arbitrary complex numbers (si,S 2 ) £ C 2 in the domain Re(s i) > —a and Re(s 2 ) > —f3. The 
corresponding Directed Polymer model was introduced and studied irf 59 . As already observed there, this model also 
contains the Strict-Weak polymer model introduced irP3 as a limit /? —> oo: 

.. ~T 7 \—n - , , P(a + rii) 

lim (/3 u)"iw" 2 = (a)„ 2 = ——-— 

>00 1 (CfcJ 

(/ 3u,v ) ~ (/3Beta(a, /3), (1 — Beta(a, /?))) ~/ 3 ->oo (Gamma(a), 1) , (35) 

which corresponds to a Strict-Weak polymer model with random Boltzmann weights on vertical edges distributed 
according to a Gamma distribution of parameter a > 0, more precisely the rescaled Boltzmann weight v! = f3u is 
distributed according to a PDF p a {u') such that 

u ~ Gamma(a) = = 7 — r{u)~ 1+a e~ u . (36) 

r («) 

A second, and completely symmetric, Strict-Weak DP limit exists for a —► 00 at fixed f3, with random Gamma{(3 ) 
weights on horizontal edges. 


3. The Inverse-Beta polymer 


We now investigate the case (ei, 62 ) = (1, —1) with /3 > 0 and a + (3 < 1, and for convenience let us introduce the 
parameter 7 as: 


7 := 1 - {a + 0) 


(37) 


In this case, a solution to the moment problem (31) is given by v = u — 1 (in agreement with the general argument 
proposed above) and u £ [l,+oo[ distributed as 


PiA u ) = 


u m v n 2 _ (_i) 


r( 7 + ff) 1 

P( 7 )P(/3) 

(Q 1 )^! (/ 3)n 2 


1 

1 - 

u 


3-1 


V = U — 1 


(a + /3) 


u £ [1, Toc[ , 7 > 0 
for m < 1 — a = 7 + f3 , ni + 712 < 1 — (a + (3) = 7 . 


(38) 


n 1 +n 2 


In this case the moments problem (311 is indeed truly solved only for n\ + ri 2 < 7 since the moments cease to exist 


beyond this bound, due to divergence for large values of u,v. However there is a more general expression of the 
moments for complex (s!, S 2 ) £ C 2 with Re(s\ + S 2 ) < 7 and Re(s 2 ) > —/3 


r 1 D d2 = 


r(7 + P) r (7 - Si - S2)T(/3 + s 2 ) 
r (7)r(/3) r(7 + ^-si) 


(39) 


Using the analytical continuation of the Gamma function to the full complex plane, one thus see, using (|39j), that 
the moment problem (31) is indeed solved in an analytical continuation sense. This situation is very similar to the 


case of the log-Gamma polymer which the present model generalizes, as we show below. Note that for the present 
model, the variable l/u is distributed according to a Beta distribution of parameters 7 and /3, 1/u ~ Beta{'y,/3) and 
for this reason we call this model the Inverse-Beta polymer. This observation renders the proof of the convergence of 
this model to the log-Gamma polymer immediate. Indeed, one has 


lim 

f3—> 00 


U V 

P'P 


2 = T( 7 - (rti +n 2 )) 

r( 7 ) 

f 1 1 — Beta( 7 , /3) 

\(3Beta("/, (3)’ f3Beta{ r ), f3) 


“'/3—>-oo 


( 1 , 1 ) 


Gamma( 7 ) 


(40) 


And this limit thus corresponds to a model of polymer with on sites Boltzmann weights (since the weights on 
neighboring links are equals in the limit) distributed according to an inverse Gamma distribution, i.e. the log-Gamma 
polymer. This analysis thus unveil a natural duality between known integrable directed polymer models, as can be 
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seen comparing (351 and (40). However, more surprisingly, this model also contains the Strict-Weak polymer model 
as a limit. Indeed, one has 


lim u ni ('yv ) n2 = 

7—>-oo 




T(/3 + n 2 ) 

r(/3) 

7(1 - Beta( 7,/3)) 


Beta( 7 ,/ 3 )’ Beta{"i,j3) 


7—>-oo 


(1, Gamma((3)) 


(41) 


which corresponds to a Strict-Weak polymer model with Boltzmann weights on horizontal edges distributed with a 
Gamma distribution of parameter /3 > 0. In terms of the partition sum, the convergence of the Inverse-Beta model to 
the log-Gamma and Strict-Weak (|2j) is easily obtained using (401 and (41). Note that one can formally take a limit on 
the moments of the Beta polymer to obtain the moments of the log-Gamma polymer. Indeed, taking on the moments 
appearing in (33) a + f3 = 1 — 7 fixed and letting a —>■ 00 , one obtain 


lim 

|a:|,|/ 3 |—>-oo,0!+/3=l—7 y —CM 


(- 1 ) 


m+n2 


(1-7) 


n 1 +n 2 


r(7 - (m + n 2 )) 

r(7) 


(42) 


But in doing so, the parameters a and /3 passes through region where the PDF Beta(a, (3) is not normalizable and 
the convergence of the Beta to the log-Gamma polymer thus does not hold in probability. The situation and relations 
between this different polymer models is summarized in Fig. [2] Notice that parts of this scheme remain empty, and 
there still remains some room for new integrable models with moments of the form (311. In Appendix [B] we attempt 
a first step in this direction by studying different analytical continuations of (31), 



FIG. 2: Duality between polymers models in the (a, /3) plane. The dashed line represents the axe of symmetry a f) (3, or 
equivalently the symmetry between vertical and horizontal edges. The blue line indicates the line a + /3 = 1 or equivalently 
7 = 1 — (a + P) = 0. Limiting polymer models are indicated by red arrows for the log-Gamma (LG) and blue arrows for the 
Strict-Weak (with weights either on horizontal edges (SWH) or vertical edges (SWV) ). We also emphasize the values of (ei, £ 2 ) 
which corresponds to the polymer considered. Notice that the region a > 0, /? > 0 and 7 < 1 is a region of coexistence of the 
Inverse-Beta and the Beta polymer, only distinguished by the value of (ei, £ 2 ). 
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IV. STUDY OF THE INVERSE-BETA POLYMER 


We now turn to the analysis of the Inverse-Beta polymer. In Section IV A we us the Bethe ansatz solvability of 
the model to obt ain formulas for the firsts (i.e. those that exist) integer moments of the partition sum of the model. 
In Section IV B we use the prescription already used in 5 -^ for the log-Gamma polymer to conjecture a formula for 
the Laplace transform of the partition sum from the knowledge of its moments. Based on this conjecture, we show 


in Section IV C the KPZ universality of the model, 
associated to the Inverse-Beta polymer. 


Finally, we study in Section IV D a zero temperature model 


A. Moments Formula and Coordinate Bethe Ansatz 


1. Coordinate Bethe Ansatz 


The moments of the Boltzmann weights of the Inverse-Beta polymer read (in the following we keep the notations 
and coordinates introduced in the general setting of the precedent section) 


u n lv n 2 = (—1) 


n 2 (oQru {P)n 2 
( OL + /3)ni+Tl2 


(43) 


where a < 0, /? > 0 and a + /3 < 1. As we showed in Section fill B[ this model is integrable using a coordinate Bethe 
ansatz Eq. (141 with a two body S'-matrix S(zi,Zj) given by (19). Its parameters are calculated from the second 


moment equation Eq. (18) and their definition (43), leading to: 


__ 1 _ 1 b_~l + a + ^_ 7 

1 + a + /3 7 —2 ’ 1 + a +/3 7 — 2 

where we have recalled the definition of the parameter 7 . We now introduce 

4 4 

7 — 1 a + f3 

Zj = e iAj , tj = itan(^) = | , zj 

Z Zj + 1 


7 = 1 — (a + P) 


1 + tj 

1 - tj - 


(44) 


(45) 


In the following we suppose c > 0, i.e. 7 > 1. As in the log-Gamma case, this is only a technical assumption that 
allows us to use the coordinate Bethe ansatz to compute the n < 7 first moments of the partition sum, and we will 
specify when the validity of some results extends to 7 < 1. Using these notations, it is a simple exercise to check that 
the S-matrix of the Inverse-Beta polymer can be expressed as 


S{zi,Zj) 


2 tj — 2 t-i + c 
2 tj -2 U-c ' 


(46) 


Remarkably, it is equal to the S'-matrix of the log-Gamma polymer studied in our previous work 5 ^. Hence, the Bethe 
eigenfunctions of this model can be taken as the one already introduced for the log-Gamma polymer, namely 


ir 


c n ) = yi Aij n 


V(a) 


A„ = 


n a 


cr£Sn 


l<a:</3 <n 


^<t(/3)) 


)• 


(47) 


Note that it only differs from the solution (20) proposed irP2 by a global multiplicative constant. Following the 
same approach than iiP^, we now study the model using periodic boundary conditions and look for eigenstates of the 
transfer matrix such that ipn( x i> ‘' 1 , Xj + L, ■ ■ ■ , x n ) = ip(xi, • • • , x n ). This imposes the Bethe equations 


e U.L = 


n 


2 ti — 2 tj T c 


i = 1 , ..n 


(48) 


Note that this is only a convenient choice and should have no effects on Z t (x) as long as t < L as discussed there. We 
will now recall some useful properties of the eigenstates (47) that were obtained^ and generalize some of them. 
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2. Recall of some properties of the eigenstates 


a. A weighted scalar product The eigenfunctions (47) form a basis of the set of periodic functions of n variables 
on N". They are orthogonal with respect to the following scalar product 


(0, i/j) = 


E 


(jeiv ,ai Tl )G{0,--- ,L— 1} ? 


IL h 




,x n )ip{x ll -■ ■ ,x n ) , h n = J| 4 _ 4 fe _ = (T- • 


k —0 


r( 7 ) 


( 49) 

b. The string solution In the large L limit, the solutions of the Bethe equations (48) organized themselves into 
strings. Each set {t a } that solve (48) is given by partitioning n into n s strings, each string containing nip particles 
where the index j = 1, • • • , n s labels the string. Inside a string, the t a are given by (we use the notations of 55 ): 


t a = tj, a = + 1 - 2a) + < E 


(50) 


k • 

where we introduced an index a = 1 , • • • , rrij that labels the rapidities inside a string, £ R denotes their common 
imaginary part and Sj iU are deviations that fall off exponentially with L. In the large L limit, the strings behave as 
independent free particles with total momentum K :] = RRl-y Ay a £ [— rrijir, rrijir). In particular, in the large L limit, 
the sum over all eigenstates can be computed as 


_ t prrij'K t poo 

E * TA [ dK :> RZ [ dk i E 


1 


1 _ + 2 
a=l 1 V 


(51) 


rrij string—states ■ J 7 * 00 

We will also need the norm of an eigenstate composed of strings in the large L limit. This was computed in 55 . 


M\ 2 = n\L n ‘ n 

l<i<j<n s 


4(fcj 

4 (ki 


7 \o _o/ \o n s I’hj rrij 

k j) + c (nij + rrij) rr mj l -pr 2 

kA 2 + c 2 (mi — m-j) 2 1 1 -f? j lJA 3 ' b 

v J n = n.= 1 h= 1 


(52) 


c. Energy-momentum of the strings Although the eigenfunctions are the same as the one for the log-Gamma 
polymer, the eigenvalues are different. The eigenvalue of the transfer matrix T n associated to an eigenstate was 
give n in |l5| as A M = RRifuf-vzf 1 ) and depends only on the first moments of the weights. Inserting their values from 
(43) and taking into account that for a string state, it is a product of string contributions, we obtain A M = IIj=i R 
with: 


a, =n< 


a 


P 1 ~ tj,a 

ol + f3 ol -|- P 1 + tj^ 

r(/3 + ^ + 2 + a^ r ( 


) = 




—/3-V-3VV + 1 


(“ 


tK o _ 7 _ Vll I i ' 
- 2 2 ' 1 ) 


ik-j 


ip- 


T( B + % 


0 r (“ 


O' 


(53) 


Where, in the second line, we have rewritten the Pochhammer symbols using Gamma functions, an indentity valid 
for integer mj. 

Another important quantity is the eigenvalue associated to the action of the unit translation operator on a string 
state, defined as 


n =n i 


i +1 


3 ,a 


- u 


r (~ir ++ i + 

r(^ + i-4) r (-^ + i + ^) 


an expression identical to the one obtained iiP2. Finally, we will also need 



2m i 


IV 

r( 


m 3 

2 


2 


i-4)r(V + | + ^) J 


(54) 


(55) 
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3. Moments formula 


We have now all the ingredients to compute the integer moments of the partition sum. As will appear clearly in 
the following, it is convenient to start with an initial condition 


Z t= o(x) — Woo5 Xt o , 


(56) 


where Wo,o is Boltzmann weight, statistically independent of the others, and distributed with an inverse Gamma 
distribution of parameter 7 . The problem with initial condition Z t= o(x) = 5 Xt 0 is obviously simply connected to this 
one and the details of the relations are given in Appendix [Cj In terms of the wave-functi on i ft(x i, • • ■ , x n ), the initial 
conditions reads ift=o{ x i, ■ ■ • ,x n ) = II”=i ^,0 an d we use the scalar product (49) to decompose it on the 

Bethe eigenstates. Using this decomposition we obtain 


i>t{x 1 , • 


o = £ 


T(y — n)n\ 

r(7)h„|| 


(A^) if^{x 1 , , Xn) . 


(57) 


In particular, 


ry , 'n. _ V t( 7 -?l)(n !) 2 t / TT ^ 

Replacing in this expression each terms by its value in the large L limit, one obtains: 


Z t (x) n = 


T (7 — n)(n\) 2 

r(7 )h n 

n s 

n^- 1 - 


n n s r +oo ii rn 3 

Ei, u 11/. ::u, 


_1_ 1 _1_ tt 4(fcj - kj) 2 + c 2 (m^ - rrij) 2 

. w-oo — n! 4 (ki-kj) 2 + c 2 (m i +mj) 2 

n s = 1 ) n .7=1 a—1 Ji a 1 <i<j<n s J J 

i 171 i 

(^TbUrT- 

1 L j,b 


-tb) 


b—1 


(58) 


(59) 


E 


rrij string— states 5 


where 


Where we have written the sum over all eigenstates as = E« fl =i 77 E( mi ... m „ ) 

E(mi ■■■ m„ )„ m eans summing over all n s -uplets (mi,-- - ,m ns ) such that EIE = and the n s \ factor avoids 
multiple counting of a same string state. Rearranging this formula and rescaling k —> ck , we finally obtain: 


Zt(x) n = n\J2^j £ A 

j=l ‘ 


n B =1 


n 771 ■ 
1=1 3 


n 


4 (hi - kj) 2 + (to,; - rrij) 2 


[ + °° dk o 

-00 27r !<i<l<Ti s “ fc j) 2 + + m l) 2 


1 ^rKE+jMM 

r(^ + i-^-) 


1 +x 


n^ + i+ikj) 


l-x-\-t 


T(/3 + ikj + ? + t=i) 


T(/3 + ikj 


_ rrij X 
2 / 


(60) 


valid for n < 7 . The convergence of the various integrals is algebraic, the integrand being 0(l/k^™ 3 ) as can be 
checked by rewriting the quotient of Gamma functions as Pochhammer symbols. This formula was checked using 
direct numerical integrations for low values of t < 2, x < 2 and n < 2. The case t = 0 is already non-trivial since it 
confirms the completeness of the eigenstates. Note that if one chooses the initial condition Z t=0 (x = 0) = <5 Xj o> then 
Z t {x) n is trivially given by (60) with an additional factor of T( 7 )/T (7 — n) in front. 


Degenerations towards the log-Gamma and Strict-Weak polymers: Since the Inverse-Beta polymer contains the 
log-Gamma polymer and the Strict-Weak polymer as limits (see (40) and (f4l])), ( |60| also contains moments formula 
for the Strict-Weak and log-Gamma cases as we show now. 


The moments of the log-Gamma polymer are obtained as the limit (Zf G (x)) n = lim^oo -X £ Z t (x) n , where 
Zjf G (x) is the partition sum of the log-Gamma polymer. Indeed, the factor exactly cancels the divergence 
of the last quotient of Gamma functions in (60), leading to the formula (54) oPE Let us recall that the present 


coordinates are t = T and i = t/2 + l as a function of those, T, X (but denoted there t, x) of that work. 


We now obtain a moment formula for the Strict-Weak polymer with initial condition Zf w (x = 0) = S x p, 

"Z t {x) n . In this case, the point-wise limit 


following 0 ' we consider the limit (Zf w (x)) n = lim. 


Ut) .... 

■7-700 r(j-n) I 
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of the integrand cannot be taken as simply and we need to first perform the change of variables kj —» kj + 
We obtain 


(Z t sw (x)) 


r r( 7 ) 

llm Tv-T 

7 -)-oo 1 (7 — n) 


7 nx n 


i y — 

h n « ! 



dkj 

~2^ 


n 

l<2<j<n s 


4 (kj 
4 {h 


A j_ m-^+7-ifc J -) V +a / r(-7 + ^) V~ z+t / r(/3 + ifc J + 7) y 
m, y r(a+ 7 -iy ; ( r(7 + ^-) j Vros + ^-^y 


kj) 2 + (rrii — rrij) 2 
kj) 2 + (m l + rrij) 2 


(61) 


Where L = — *3) + M. Since the integral over kj quickly converges as O^l/k 2 ™ 3 ), we can now close the different 
contours of integrations on the upper half plane before taking the limit 7 -7 00 . This leads to: 


{Zf^{x)Y=n\ y 

n s = 1 


1 

n s \ 


E 


(mi,.. 


n 

)» • 


dkj 

2tt 


n 

l<i<ji<n s 


4(fcj 

4(fci 


ft 7 /rX-7+^)y-" +t / r(/3 + T% + 7) V 

jimj \ T^+ih) ) \T(P + ikj-?t)) 


kj) 2 + {rrii - rrij) 2 
kj) 2 + (rrii + mj) 2 


(62) 


where L is an horizontal line that stays below all the poles of the integrand. This formula is formal because the 
resulting integral does not converge, but one must remember that we have formally already closed the contours 
of integrations. Computing the integral on ki thus just amounts at taking the sum over the residues of all the 
poles of the integrands except those of the type ki = kj — iA where A > 0 (since the contours have been closed 
on the upper half-plane). 


B. Fredholm determinant formulas and KPZ universality 


In this section, we use the formula (60) to obtain the Laplace transform of the distribution of Z t (x ), 


gt,x(u) =exp (-uZ t (x)) . 


(63) 


The issue of obtaining this generating function from the sole knowledge of the integer moments of the partition sum 
was thoroughly discussed in 22 and here we follow the same route. 


1. The moment generating function 


We start by computing the moment generating function 

97; 


(“) = £ 


(64) 


n —0 


where u > 0. Here, though Z t (x) n is only defined for n < 7 , the right hand side of formula (60) is well defined for 
n £ N (except if 7 £ N) and we take advantage of this analytical continuation to perform the sum (64). Note that 
this object has no reason to correspond to the Laplace transform of Z t (x) but 67 , as in the log-Gamma case, we will 


use it to conjecture a formula for the true Laplace transform g t x (u) defined in (63). Since we perform the sum over 


n £ N, the constrained sum appearing in (60) becomes free summation and one can write 


mom 
Vt 


+00 1 

( 11 ) = 1 + y — -Z(n s ,u) 

L ^ n.J 


(65) 


67 And indeed it is not, a simple reason being that, just as the Laplace transform of the PDF p 7 ,/3 of the Boltzmann weights of the 
Inverse-Beta polymer (see |38j), the Laplace transform of Zt(x) is not an analytic function. See also Appendix |d| for more details on 
this question. 
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where 


n s +oo -+00 77, 

z(n s , U ) = n e / of n 


i=i m j—1 ■ 




4(fcj — kj) 2 + ( nrii — rrij ) 2 
4(fcj - kj) 2 + (rrii + mj ) 2 


n 

1=1 




1+x 


r(^ + i + ifc,) 


1—:r+£ 


r(/3 + ikj + \ + si) 
r(/3 + ikj + 2 _ 


In this formula and following 5 ' 5 !, 0 ne recognizes the structure of a Fredholm determinant 

5 ™°”» = Det + 


with the kernel: 


Ts-mom 
*^t,x 


_°°_ /*+°° fib. 

(v 1)W2 )=E / _(- u ) m e - 2ife(,,1 - , ' 2) - m(, ' 1+U2) 

m= lJ-oo TT 

Uos + ifc + i-f) 


r (-7r + J-^)V +x /r(-f+ | + *fc)V x+t ^r(/3 + ifc + 2 + ^) 




r(f + |+*fc) 


( 66 ) 


(67) 


( 68 ) 


and A'™ om : L 2 { 


_), so that the two auxiliary integration variables "17 and V 2 are positive 68 . 


2. The Laplace transform as a Fredholm determinant 


We now use the same prescription used in 5 -® to obtain a conjecture for the Laplace transform gt,x(u) from the 
moment generating function g”“ m (ii). It consists in rewriting the sum over m in the Kernel K™° m as a Mcllin-Barnes 
integral. In Appendix [D] we also show how this type of manipulation efficiently works on a simpler object, namely 
the Laplace transform of the PDF p 7)/ 3 defined in (38). We thus conjecture, gtx{u) = Det (/ + K tx ) with 


Kt,x(v i,i> 2 ) = 


r ’ +tx> dk -1 

7T 2 i 


ds 


c sin(7rs) 

r C-| + l-ik)^ 1+x 
r(f + %-ik) 


u s e -2ik(v 1 -v 2 )-s{vi+V2) 


r (-f + i 


■ ik) 




ik) 


1 -x-\-t 


r(/3 + ifc + ^ + f) 
r{p + ik + %- f) 


(69) 


where (7 = 0 + ® with 0 < a < min(l, 7 ) and A7, x : L 2 (R + ) -+ A 2 (R + ). As in the log-Gamma case, we expect 
this formula to be also valid for 0 < 7 < 1. Note that in going from ( 68 ) to (69) we have to choose an analytical 


continuation to go from m € N to s £ C. Here the chosen analytical continuation is the most natural one in the sense 
that it generalizes the one used for the log-Gamma polymer iif®, and also mimics the calculation of Appendix [d] 
This Kernel is the one that is naturally obtained from the Bethe Ansatz and its structure is reminiscent of the string 
solution: the integral over s encodes for the contributions of the different types of strings, whereas the integral over 
k is the summation on the momenta of the strings. As shown in 1 ® (section 11 ), it is also possible to rewrite gtx(u) as 
the Fredholm determinant of another Kernel which contains one less integral. Since the proof is strictly analogous to 
the case of the log-Gamma polymer, we only give here the final result: we also have gt, x {u) = Det (/ + K^ X A ) where 


K^{z,z')= [ 

J 2, 


dw 


2a+a+* 47 t(w - z') sin(7r(w; - z)) 

r(y + a — w)\ 1+x f T(z — a) \ 1 x+t f T(w — a + /3)' 4 


r (7 + a - z) 


F(w — a) 


T{z - a + 0) 


(70) 


where : 

should be valid for arbitrary x, whereas for the log-Gamma polymer the analogous formula was only valid for 2 x < t 


L 2 (a + a + *1 


L 2 (a + a + ®) 0 < a < min(l, 7 ) and 0 < a < 7 — a. Note that here this formula 


68 Note that this FD structure would have been broken by the initial condition Zt=o(x ) = S x= o■ 
factorizable term of the form r( 7 )/r (7 — m «) )• 


(In which case ( |66| ) contains a non 
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(with a mirror formula for the other case) 69 . This formula is a large contour formula and an analogous small contour 
formula should also exist, as in the log-Gamma polymer. Let us also mention here that, following the sam e p roce dure 
that led in the log-Gamma case to formula (63) and (64) of^, it is possible to directly obtain from (69) or (70) 
formulas for the PDF of log Z t (x) as differences of two Fredholm determinants. 


3. The Laplace transform as a n-fold integral 


IiPl, a formula giving an identity between a certain class of Fredholm determinant with Kernels similar to the one 
in (70) and a class of n-fold contour integrals was given (Theorem 2). Though the explicit form of (70) explicitly 
breaks the hypothesis under which this formula was proven, an analogous formula should also exist in a more general 
setting. Guided by this belief, we conjecture the following formula for the Laplace transform: 



with 0 < a < min(l, 7 ), valid for Re(u) > 0, 1 < J < I and where x = / — 1 and t = I + J — 2. This can be seen 
as a modification to our model of the formula given irt®2 (Theorem 3.8), and also stated irP^ (Proposition 1.4), for 
the log-Gamma polymer. Since ours is merely a conjecture, we have tested it numerically against direct numerical 
computations of the Laplace transform for various u, (3 and 7 and for J = 1 , I = 1, 2, 3 and J = 2, I = 2. 


4 . Degeneration towards the log-Gamma polymer. 


The results of the last three paragraphs for the Laplace transforms of Z are easily seen to degenerate into the 
usual results for the log-Gamma polymer as (3 —► +00 using that exp (—uZ^ G (x)) = lirng^oo exp(— jcZ t (x)). For 
example, taking the limit on formula (69), this introduces a term exp(— st log(/3)) in the Kernel that exactly cancels 
the divergence of the last quotient of Gamma functions, and similarly for the other formula. 


C. The large length limit and the KPZ universality. 


We now study the limit of polymers of large length 1 > 1 for polymers with fixed endpoints (0,0) and (t, x) = 
(f, ( 1/2 + if)t) where ip £ [— 1 / 2 , 1 / 2 ] represents the average angle of the path measured from the diagonal of the 
square lattice. The large t behavior of (69) is estimated through a saddle-point analysis similar to the one irP^ to 
which we refer for details. We define 


,1 


G<p(y) = G + ^)!ogr(^ -y) ~G~ ^)iogr(^ + y ) + io g r(/3 + 7 + y ). 


,1 


7 


7 


(72) 


So that the leading behavior of the product of Gamma functions appearing in (691 is 


(Fr)* := exp (t (G v { S -+ik) - G v {- S - + ik))) . 


We now look for the critical point (s, k) = (0, — ik v ) such that G”(k v ) is 0. This defines implicitly k v as 


(73) 


{\ + vW'Q - M + M + 1 + M = o . 


( 74 ) 


69 Convergence of the w integral is checked using that \T(x + iy)\ — y/2n\y\ x 2 e 2 1^1. 
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Where if = U- is the diGamma function . Expanding (73) around this critical point, one obtain 


(IT)* = exp ( t ( G'Jk v )s + 


<%'(**) S q„7.2n 


6 


(-3 sk 

4 


Q(s 4 


(75) 


where k = k + ik v and s are considered to be of the same order (this is consistent with the rest of the calculation, see 
below). The linear term G'(k v ) corresponds to an additive constant in the limiting free energy, whereas the cubic 
term sets the scale of the free-energy fluctuations. To pursue the asymptotic analysis, we define 

F t (ip) = - log Z t (x = (1/2 + tp)t) = c v t + X V ft{v) 

'tG"{k v ) 


c v = -G’^kv) 


A v — 


8 


9 t,<p(z) = exp (—e -’vG+AM)) (76) 

Where F t (<p) is the free-energy of the directed polymer and gt,<p(z) is a rescaled Laplace transform which has a proper 
t —> oo limit for fixed zel. Indeed, since gt,x=(i/2+ip)t{ u ) can be written exp (—e log (") _Ft ^)), one has the identity 
gt v {z) = g t x =(i/ 2 +<p)t( u = e Cvt ~ XvZ ). Rescaling s — > s/X v , k —> Vi —> A v Vi and inserting u = e c v t ~K z j a s well as 

0 

_ . r , > , ^ . r o „3 

(77) 


the expansion (75), into (69), one obtains gt, v (z) = Det + K tttp (v 1 , 1 V ' 70 

ds 


~ f dk -1 

Ft,<p(Vl,V2) = / —— 


— sz — 2 ik(v 1 — v 2 ) — s(vi+V2)— 4 fe 2 s+^-+ 0 ( W) 


7r 2 i J c X v sin(7rA-) 


where K t>lp : L 2 (R + ) 


L 2 (R + ) . The large polymer length limit X v —» oo can be safely taken in this last expression, 
leading to a kernel for which there is more freedom in the choice of the integration contour C: it should only 
define a convergent integral and passes to the right of zero. The f —» oo limit of the rescaled generating function can 
thus be written as linp^oo gt, v (z ) = Prob(—f < z) = Det(I + K 0 0 ) where I\ ^ : L 2 (R + ) —> L 2 (R + ) is given by 


Koo(vi,V 2 ) = - f ^ 


dyAi(y + z + v\ + v 2 + k 2 )e l ^ Vl ^ 


(78) 


where we used the Airy trick J R dyAi(y)e ys = e^~ valid for Re(s) > 0, followed by the shift y —¥ y + z + v\+v 2 + 4 k 2 , 
the identity f c e sy = 6(y ), and the rescaling k —t k/2 . As irP^, this kernel corresponds to the Tracy-Widom GUE 
distribution as det(I + K 00 ) = F 2 (2~^z) where F 2 (z) is the standard GUE Tracy-Widom cumulative distribution 
function. We have thus shown 


lim Prob 

t—> OO 


logZ t ((l /2 + y)f) +tc v 

X in 


<232 ) = F 2 (z) 


where the ((^-dependent) constants are determined by the system of equations: 

o = - K) ~ i\ ~ + M + ^'(Z 3 + \ + K) 

c<p = (~ + - M + (^ - + k v ) - tH/3 + \ + M 


t 


1 


K = ( ((h + vWd, - M + (^ - ¥#"(77 + K) - i + K) 


i 


r'y 


(79) 

(80) 
(81) 

(82) 


Angle of maximal probability The free energy per unit length c v is maximal in the direction defined by the angle 
< p * such that = 0. It is easily seen from (80) that it is realized for k v = 0, and ip* is thus given by 

1 V> ; (/3 + 7/2) 

'2 ^(7/2) 


y?* = 


<0 , 


(83) 


70 


The extra factor e 2k <p x <p( v i v i ) originating from the change of variable has been removed since it is immaterial in the calculation of 
the Fredholm determinant. 
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and the optimal energy per unit length is thus 

c* = c v * = V>(7/2) - i>(P + 7 / 2 ) • (84) 

The amplitude of the fluctuations in the direction ip* are 




1/3 

= (j^'(0 + 7/2)-<(7/2)) 

—/3^o Qv ,, "(7/2)/3^ 

—/3—voo (-^"(7/2)) • 

And one recognizes the usual log-Gamma result for <p = 0. In the log-Gamma limit /? 
the parameter (3 > 0 biases the DP towards the vertical direction. More precisely, 


(85) 

00 , one recovers p* = 0 , but 


ip* ~/3—>0 


¥>* —/3->oo 


^(7/2) 

^(7/2) 


1 


2^( 7 /2)/3) 


£ + 0 (/I 2 ) 
0 ( 1 /^ 2 ) . 


( 86 ) 


For small displacement around this optimum direction ip = <p * +5p, one retrieves an i sotr opic continuum limit 
characterized by an elastic coefficient tt such that c v ~ c^,* — jn5p 2 . One easily find using (801: 


k = —t 


W>'(7/2)) 2 


^"( 7 / 2 ) - ip"(/3 
which generalizes the known result for the log-Gamma. 

Degeneration towards the log-Gamma and Strict-Weak polymers. 


7/2) 


(87) 


The Laplace transform of the partition sum of the log-Gamma polymer is obtained as exp(— uZpf G (x)) = 
lim^oo exp(— jiZ t (x)). This amounts to change u s —> u s exp(— ts log(/3)) in the above formulas. For large (3 we 
use the limits ip'(x) —> x ^oo 0, ip”(x) — > x ^oo 0 and ip(x) = log(a;) — ^ + O(^). It is then easily seen that the 
presence of f3 do not change the position of k v in this limit nor the amplitudes of the fluctuations A^,, whereas c v 
receives a contribution proportional to — log(/3) which exactly cancels the rescaling of the partition sum. This 
shows that the system of equation (801 converges to the one of the log-Gamma. 


• In the case of the Strict-Weak polymer, the rescaling of the partition sum introduces a term that amounts 
to change u s —» u s exp(si log(y)) = exp(s(l/2 + y>)log( 7 )) in the above formulas. This suggest to look for a 
solution of the form ^ + k^ w . The system of equation (80) then converges to 


0 = -b-^ , (O + V' , (/3 + O 


4 W = U ) - W + K w ) 


,1 


SW\ 


'„SW\ 


kSW 


= \-i\(n-'pW(K w )-r(P + K w ) 


( 88 ) 

(89) 

(90) 


so that we retrieve the result oiP^ for the Strict-Weak polymer case (the precise correspondence with their 


notations reads k = l/(l /2 -ip),t = k v , k = [3, f k , K = ~nc^ w and g k = t(1/ *_ v) (A y )- 


D. A low temperature limit. 

1. Definition of the zero temperature model 

In this section we study the limit 7 = ey' and /3 = ef3' of the model with e —> 0 (hence, a —i 1). As we show now, 
this model converges to a zero temperature problem. 
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The analysis is similar tcP^. There (Lemma 4.1) is was shown that for a random variable z chosen with a Beta(a = 
ea , 0 = eb) distribution, the joint PDF of the pair (— elnz, —eln(l — z)) converges in law to (£E a , (1 — £)Eb) as e —> 0 
where £ a Bernouilli random variable (i.e. £ = 0,1 with probabilities p = b/{a + b), 1 — p) and E a ,Eb exponential 
random variables of parameters a and b respectively (i.e. p(E) = ae~ aE 6(E)) statistically independent from £. Note 
that the correlations between E a and Eb are unimportant since they are multiplied by £ and 1 — £ which cannot be 
non-zero simultaneously. The occurence of the Bernouilli variable is intuitively understood since in that limit p{u) 
exhibits two peaks, one near u = 0 and one near u = 1 with weights p and 1 — p, and the exponential distributions 
arise by zooming-in on these peaks and rescaling (u for the first peak, v = 1 — u for the other peak). 

Since in the Inverse-Beta model \/u is distributed as a Beta("f, 0) random variable, we immediately obtain that 
the rescaled random energies of the model (£„,£„) = (—elog(u), —elog(u)) converge in probability to 

(—elog(u),—elog(u)) (~(E y , (1. - £)£> - (Ey) = ( 91 ) 

where £ is a Bernoulli random variable of parameter p = 0 , /("f r + 0'), Ey and Ep> are exponential random variables 
of parameter 7 ' > 0 and 0' > 0, independent of £. Equivalently one can choose: 

(£(,,£') = (0, Ep’) , with proba 1— p (92) 

{£' u ,£' v ) =—Ey{ 1,1) , with proba p (93) 

i.e. a model where disorder is chosen randomly either on the site or on the pair of edges arriving at it, with a penalty 
for the horizontal edge. The two cases corresponds to two peaks near (u,v) = (1,0) and (u,v) = (+ 00 ,+ 00 ) in their 
distribution in that limit. In terms of the partition sum of the polymer, the limit reads 

- elog(Z t (a:)) =-elog( Y exp(J^ log(w e ))) 

7r:(0,0)->(t,x) ee-7T 

= —elog( 

7r:(0,0)->(t,x) ee-ir 

~e_>o min T:(0;0 ^ (t;X ) Y, K. ■= £(t,x) ■ (94) 

eE7r 


which justifies the name zero temperature limit: the rescaled free energy of the original model converges in probability 
to the minimal energy £(t,x) f° r the set of all polymers with starting points ( 0 , 0 ) and ending points (t, x ) in the random 
environment with energies £' e distributed according to (91). 

Degeneration to the Exponential (i.e. q = 1) Johansson model: In the so-called log-Gamma limit, i.e. 0' —> + 00 , 
one obtains p = 1 hence: 


(£' U ,£' V ) = -E y ( 1,1) (95) 

i.e. the on-site exponential distribution model of parameter 7 , also identical to the q —> 1 limit of the Johansson 
model, studied iiM Note that the extra weight wqq at the origin which we included, allows to precisely recover 
the Johansson polymer model (with an exponential variable also on the site x = t = 0.). To make contact with the 
notations of 21 I we have = —H(M, N) with M = I = 1 + x and N = J = l+ t — x. 


Degeneration to the zero-temperature limit of the Strict- Weak model: In the limit 7 ' —> 00 one obtains p = 0 hence: 

(CO = (0 ,Ep) (96) 

This model can be interpreted as a discretization of a zero temperature version of the semi-discrete polymer model 
where one replaces the set of independent Brownian motions by a set of independent random walks. 


2. Fredholm determinant formula for the zero temperature model 

In order to obtain a Fredholm determinant for the zero temperature model starting from our expressions for 
g tx (u) = exp(— uZ t (x)), we rescale u as u = exp(r/e) with r G R. fixed. Indeed, one then has 

gtx{~ exp(r/e)) = exp(-exp -(r + e\ogZ t (x))) 

e 


—■ ►e->o 0(—r + <£(t, x )) 
= Pro6(€ (t) x) > r) . 


( 97 ) 
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We can thus directly write a Fredholm determinant formula for Prob((£(t,x) > r ) by inserting u = exp(r/e) in (69l. 
The point-wise e —» 0 limit of the Kernel is taken using a rescaling s —> es, a — > ea (so that the contour C do not 
crosses poles when we take the limit), k —> ek and v z —i v^/e and using 


—y ~~K-s-o — > r(ea;) — £ _>.o-k 0(1) . 

sm( 7 res) ns ex 


We thus obtain Prob> r ) = Det (/ + Kj~°) with 

/ ' +0 ° dk [ ds 


K t,x( v 1 W 2 ) = — 


sr—2ik{y\—V2)—s{v\-\-V2) 


n r 2 ins 


(98) 


(99) 


/ \ l-\-X / t \ 1 — X-\-t / / \ t 

\ — ik \ ( § + \ + ik \ (/3' + ik — § \ 




-f + £ + ** 


\P' + + "2—k§y 

where now C = a + with 0 < a < 7 ' and Kf x ° : L 2 (K + ) —>■ L 2 (R + ). Using the same type of rescaling as above, 


we also obtain an analogous expression to (70) as Prob(<£( t ,x) > r ) = Det (^1 + K^ X A,T where 


K, 


BA,T -0 
t.x 


M= [ 

j 2 


dw- 


p r(w-z) 


2a+a+m 4n(w - z') n{w - z) 


/ 1 \ 1+* / \ 1 —x+t / | Ol \ t 

7 + a — z \ / w — a \ ( z — a + p 


7 '+ a — w 


w — a + ft' 


( 100 ) 


where : K^ X ’ T 0 : L 2 (a + a + «K) —» L 2 (a + a + iK) 0 < a < 7 ' and 0 < a < 7 ' — a. 

We also immediately obtain a formula analogous to our conjecture as a conjecture for the T = 0 model: for 

I > J 


1 f A duu A A e r ^- Q ) 


Pr< * (<E( ‘'- |>r)= S] W ,P 1 ]J., |,r " l ,!!(-»,)' VY + «-»7 V^-a + Y 


7' 


7+J-2 


• ( 101 ) 


with 0 < a < 7 '. 

Limit to the Johansson model: For ft'/"/' = +00 one thus finds a formula for the DP model of Johansson (i.e. with 
independent exponentially distributed on-site energies). It is then interesting to compare our formula with the one 
obtained irPH (formula (1.18), which reads (for r < 0), I > J > 1: 


N 


Prob(—£( t , x) < ~r) = -gr / Yl dx o II ( x * “ x i) 2 II 


IJ J[0-rY j=l 1 <i<j<J 


x I ~ J e~ x i 


( 102 ) 


l=i 


which coincides with the CDF of the largest eigenvalue of the Laguerre Unitary Ensemble (LUE) of random matrices 
(the constant Z\j simply ensures the normalisation to unity of the measure on (R + ) J ). 


3. Asymptotic analysis and KPZ universality for the zero temperature model 


We now study the large length limit of the zero temperature model: t —>■ 00 and x = (1/2 + ip)t- The analysis 
is similar to the one made for the finite temperature model and here we only give the main steps. As before, the 
t -> 00 limit is dominated by a saddle point. The dominating term in the Fredholm determinant (99) now reads 
exp(t(G v ,(| + ik) - G v (-§ + ik))) with 


,1 


G*{y) = + v) lo g(-7r ^ v) + (a - <fi) lo §(~o" +y)~ iog(/3' + /r + y) ■ 


,1 


7' 


7' 


(103) 


Note that with the contour previously chosen the arguments of G^ stay away from the branch cut of the logarithm. 
As before we look for a critical point, (s, k) = (0, —ik v ) such that G’^(k v ) = 0. This defines k v as 


(§ + f) 


(h~<p) 


1 


(1-M 2 (1 + M 2 C 8' + l + k v y 


= 0 . 


(104) 
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Note that this equation as in general several solutions, but the only physical one must have \k v \ < y'/2 to truly 
dominate the integration. To this point, we can now follow the exact same steps as before by taking 


r = tc lp - X v z 


c-ip G v (ie v ) , A 


K'Ck) 


(105) 


and using the same rescalings in (99). In the large length limit, this leads to 

{l/2+tp)t) tc<p 


lim Prob 

£—>■00 




with 


C(£> - 


(h + P) (h-v) 


> -2a z = F 2 {z ) 


1 


“ 2 — fcw "t" + kp [F + -)*- + k, 


0 = 


2 2 1 n( P 

(s + v) i\-v) 


2 1 "'V 

1 


(t ^ M 2 (t + M 2 (/?' +1 + M 2 


* = / t I (l + 2 y) | ( 1 - 2 ^) 


(i-M 3 (y + M 3 (/3' + y + M 3 


(106) 

(107) 

(108) 

(109) 


Note that this result is coherent with the one obtained at finite temperature (79) and (80) and can be obtained from 
it by scaling 7 = ey', /3 = e/3' and k v = ek v . 

Angle of optimal energy The angle of minimum energy ip* of the model is obtained by solving = 0. This 

imposes k v = 0 and, using (107), we thus obtain 


y 2 1 

05 * =--- 7 - < 0. 

8 (y + ^)2 


( 110 ) 


As for the finite temperature model, we thus retrieve that /3' > 0 biases the DP towards the vertical direction. For 
/3' —»• 0 we obtain once again <p* = — A. The optimal energy per unit length, and the scaling parameter A v * at the 
optimal angle are respectively 




2/3' 


A( 0 * — 


YU3' + t) ’ 


(y) 3 (y - 1 - 2/3') 3 


1 3 


(111) 


recovering the results for the Johansson model: 


In the limit /?' = +00 the above equations (107) can be solved explicitly. One finds k v = — ^(1 — \/l — 4</> 2 ), 


where we have chosen the root which vanished at the optimal angle p* = 0 (i.e. the diagonal which is a symmetry 
axis in this case). This yields: 


—-/(l + \/l 4<p 2 ) , \<p — — t 3 


8 p 4 


y ' ■ y \(i-yi-4 <y)Vi-4<y 

We can now compare with Johansson result (formula 1.22 iiP3) which reads (for y' = 1 ): 

H(gJ, J) ~j_> +00 (1 + Jg?J + <r 1/6 (i + VaY /3 J 1/3 X 2 


( 112 ) 


(113) 


where \2 is a Tracy-Widom GUE random variable (of CDF given by F 2 ). With a little bit of algebra one can check 
that this is exactly equivalent to our result, namely: 


^(£,x=(l/2+y?)£) —£—>•+00 2 / A, 

with £(t,x=(i/ 2 +ip)t) — -H(gJ, J), taking into account that J = 1 -\-1 — x ~ hence g = 


(£,x=(l/2+<p)£) —£->--|-oo Z ' AtpX 2 (H4) 

1 XW - 1+2 7 
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V. CONCLUSION 


In this paper we attempted a classification of finite temperature directed polymer models on the square lattice with 
homogeneously distributed random Boltzmann weights and a certain type of short-range correlations (Section IIAI, 
for which the moments of the partition sum Z t (x) can b e c alculated via a coordinate Bethe ansatz. Following the 
pioneering work oiP^l we obtained a rigorous expression (29) that constrains the possible forms for the moments of 


the underlying distribution of weights. We discussed in details the possibilities of finding PDF’s with the appropriate 
moments (29) and, though the classification is still not complete, we were able to exclude a large number of cases. In 


cases where the moment problem has a solution, we retrieved all the previously known finite temperature integrable 
DP models (Section III B), and introduced a new one, the Inverse-Beta polymer, which appears as a natural two 
parameters generalization of the log-Gamma polymer, but also contains the Strict-Weak polymer as a limit. Using 
the Bethe ansatz, we obtained an integral formula for the moments of the partition sum (60) of the Inverse-Beta 
polymer, with point-to-point boundary conditions. Along this route, most of the tools developed irP^ for the Bethe 
ansatz solution of the log-Gamma polymer proved very useful and were generalized. 

Starting from the moments formula and using analytical continuations, we obtained two equivalent Fredholm 
determinant formulas for the Laplace transform of the PDF of the partition sum (69) and ( |70| ), and conjectured 
a n-fold integral formula for the same object, which generalizes a known formula for the log-Gamma polymer 
obtained in 53 in the framework of the gRSK correspondence. Using our Fredholm determinant formulas and an 
asymptotic analysis in the limit of large polymer length, we were able to obtain the KPZ universality of the model 
(critical exponents and Tracy-Widom GUE free-energy fluctuations) (79) and as well as exact implicit expressions for 
the mean free energy and the amplitude of fluctuations as a function of the polymer orientation w.r.t. the diagonal. As 
an application we obtained an exact expression for the optimal angle which minimizes the free-energy of the polymer 

PI- 


In Section |IVD| we introduced a zero-temperature DP model as a limit of the Inverse-Beta polymer, which generalizes 
the previously known zero-temperature limit of the log-Gamma polymer. Using the exact formulas obtained for the 
Inverse-Beta polymer, we showed analogous formulas for this zero-temperature model. In particular we obtained 
exact formulas (Fredholm determinant and n-fold integrals) for the cumulative distribution of optimal energy of the 
model (99), ( 100|) and (|101[). Using an asymptotic analysis, we showed the KPZ universality (106) of the model. Our 


formula compare successfully with some results obtained by Johansson in his pioneering study of the Exponential 
zero-temperature polymer^, a particular case of our zero-temperature model. 

We believe that the present work could be used as a guide for future research of new integrable DP models. In 
addition, we once again showed that the replica Bethe ansatz method is a valuable and versatile tool for the analysis 
of such DP models. In particular, some results of this paper could prove useful and adaptable to the analysis of the 
model with different boundary conditions. 

For future works on the Inverse-Beta polymer, it should be very interesting to obtain a solution of thi s mo del using 
the gRSK correspondence or a generalization of the latter (as in the recent worlP2l). Our conjecture ( |7lj ) could be 
proven (or invalided) using these techniques. In addition, we know that an inhomogeneous version of the log-Gamma 
polymer was amenable to analytical treatment in the framework of the gRSK correspondence, and it is thus likely 
that an inhomogeneous version of the Inverse-Beta model should also exist. 

For future works on the classification of directed polymer models, various directions of research remain. The most 
direct one is to understand if some integrable models remain to be found to fill the left voids in Fig. [2] (as e.g. our 
proposal of Appendix [B]). Other directions would be to extend this framework to introduce inhomogeneous models, 
or different disorder correlations. The precise implications of our classification of finite temperature DP model for 
possible zero temperature integrable DP models remain to be elucidated. Indeed all the models we found in our 
framework admit a zero temperature limit. For example, the zero temperature limit of the log-Gamma model is the 
q —> 1 limit (beware that this q is a priori different from the one used in Section III) of the zero temperature model 
of Johansson^, i.e. the Exponential zero-T model (as was pointed out iiP). However, at this stage, our framework 
seems to miss the q ^ 1 case of the Johansson model. A natural question is then to understand if a finite temperature 
integrable DP model sits above the Johansson model Vg, and whether the zero temperature model studied in this 
paper admits a g / 1 generalization. Since Johansson’s model is determinantal, a related outstanding question is 
to obtain a deeper and more systematic understanding of the relations between Bethe ansatz solvable models and 
determinantal processes which seem to often occur as limit cases of the former. 
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Appendix A: The |g| < 1 case: study of degenerations. 


Here we study in details the possible degenerations of the parameters (q, v, p) that would eventually lead to a PDF 
p(u,v ) such that the moments of u and v are given by (29) and n max > 2. As in the main text, we restrict to the 
domain |g| < 1 and consider the random variable z x = u + xv, x € R. Its variance is: 


-~ 2 C = (p- 1)(1 - q)(p ~ v){px - 1 ){px - v) 


p 2 {v — 1 ) 2 (vq — 1) 


(Al) 


which must be positive. Since the polynomial in x changes sign at x = ^ and x = f one must look for cases where 
/ = j l ■ The different cases to investigate are thus v = 1, |/z| = oo, p = 0 and a combination of these cases. It is 
instructive to look at the variance for x = 1: 


2 C (m- 1 ) 2 ( 1 - q)(p- 


zt = 


p 2 (v — 1 ) 2 (vq — 1) 


(A2) 


One sees that the positivity of the variance implies qv > 1, as long as we do not consider degenerations p 
or q —> 1 (in which case the variance of Z\ may vanish and the condition may disappear). 


1, p 


• If v -A 1, it is easy to see from the variance of z\ that there must be at least one additional degeneration, 
either (i) p —> 1 or (ii) q —> 1 (iii) both q —»■ 1 and p —> 1. The first one can be ruled out as follows: setting 
p = 1 + ae and v = 1 + be with e —> 0, one finds that the variance is z 2 ° = a{ - b ~ a > (1 _ x) 2 and that from (29) the 


marginals have integer moments u n = 1 — f and v n = | for all n > 1. This implies that u + v = 1 and v = 0,1 
with probability #, which then predicts joints moments different from the ones obtained from (29) in that limit, 
hence no joint PDF exists in case (i). The case (ii) and (iii) both imply a q —> 1 limit which we discuss in the 
end of this appendix. 

• If \p\ -A oo. In this case, looking at the original moments ( [29] ) , we see that we must scale v -A v/p to obtain a 
well defined random variable. We define (u',v') = ( u,v/p ). In the limit p -A oo, the moments are 


< j f7~L\ qjf'Tl'2 — 


(-1 ) n2 q^s 


2 O 2 - 1 ) 


( 9 ; q) 


n 1 +n 2 


1 


iy\ < l)n 1 +n 2 (q;q)n 1 (q- 1 q) rt2 ^n\-\-n 2 

we now define the random variable z' x = u' + xv 1 and compute its variance: 

(1 - q)(x - l)(x - u) 


( 4) 2 = 


{v-l) 2 {vq - 1 ) 


(A3) 


(A4) 


which must also be positive for all x. However, this polynomial changes sign at x = 1 and x = v so we must 
have v = 1. Since the constraint qv > 1 still holds and q < 1, the only possibility is to have q A 1 as well, a 
case discussed below. 


If p = 0, looking at (29), we see that we must now rescale (u,v) as (u',v') = (pu,v ) and we obtain 

n l(ni-l) 


y/ n l'l)? n 2 = 


(-1 ) ni q 


( 9 ; q) 


n 1 +n 2 


1 


(y\ q)n 1 +n 2 ( 9 ; 9)711(9; 9 ) 712 ^ 711+772 

i.e. this is a case identical to the previous one, and we also conclude that we must have q —> 1 . 


(A5) 


Let us now discuss all the possibilities in the q —> 1 limit. Taking the limit directly on (29), one obtains 

(l-v/p) ni (l- p) 712 


u n\ v n 2 _ 


(l-*) 


771+772 


(A6) 


where we used that at fixed n , a, ( q a ; q) n (1 — q) n {a) n , where ( a) n = a(a + l)..(a + n — 1), and here we took 

a = 1. Obviously, the limit we took only works if u, p and v/p are all different from 1, but it also encompasses 
other limits such that the p —» 00 and q —> 1 case discussed when analyzing (|A3|) . The moments (A6) correspond to 


deterministic weights u = 1 and v = These models are obviously integrable, but trivial. 

We must thus study the q —> 1 limit with at least one of those parameters that goes to 1. The question of the speed 
of the convergence then arises. In general, taking q = 1 — e and a = 1 — a'E with £ > 0, one has (a; q) n — g ->o £ n (a')n 























26 


(if £ = 1), ( a\q) n — q ^o e ? a'e n_1 (n— 1)! (if C > 1 and n > 1) and ( a;q) n e n ^{a') n (if £ < 1). The possibility 

of using £ < 1 for the convergence of p and/or v (slow convergence compared to q) is uninteresting since it leads 
to pure power-laws. The possibility of using £ > 1 (fast convergence compared to q) is also uninteresting since one 
cannot rescale (u,v) to obtain well-defined moments in the e —> 0 limit. In the following, we thus only consider the 
possibility of convergence of the parameters at the same speed than q. Let us first examine the cases where only one 
of those parameters goes to 1. We obtain 

• If p = qP and v 1, 


u Ul v n2 = (/3)„ 2 x power — laws 


(A7) 


If v = q a+ P and p^ l, 


U n i v n 2 = 


(ck + /3)m+n 2 


x powei — laws 


(A8) 


If v/p = q a , v 7^ 1 and /i / 1, 


u n lv n 2 = (a) ni x power — laws 


(A9) 


where we have not written the precise form of the unimportant power-law terms. As discussed in the main text (Section 
III B), these cases indeed correspond to proper PDF’s and known integrable models for some range of parameters a, 
f3. The first and third ones correspond to the moments of the Strict-Weak polymers, the second one corresponds to 
the moments of the log-Gamma polymer. Notice however that these models can all be obtained as limits of the case 
(g, p , v) —> (1,1,1). Indeed, taking p = q 13 , u = q a+ P and q —> 1, we obtain 


u n ly n 2 = 


(oQni (fi)n 2 

(o + f3) ni -\-n 2 


(A10) 


Taking ap propr iate limits on this last formula, we can retrieve the three precedent cases. For example the limit 
/3 —> oo of (A10) leads to, rescaling u as f3u, ( (3u) ni v n2 = (a) ni , and we obtain the moments of the third case. Taking 
a similar limit |a|, |/3| oo with a + (3 fixed, we obtain the second case. 


Hence, in this sense, the most general limit of (291 that can lead to distributions with a well-defined variance is the 
(g, p , v) (1,1,1) limit. This limit is studied in details in Section IIIB 


Appendix B: A more systematic study of analytical continuations of moments. 

Starting with a polymer model with moments given by 


u n lv n 2 = ( ei )"i ( e2 ) r 


(ck)rai (/3)n 2 

(ck + /3) ni+ „ 2 


(Bl) 


where (ck,/3) are arbitrary, (ei,e2) € {—1, l} 2 , we look for distributions p{upv) such that (Bl) corresponds to the 


moments of positive random variables. It is natural, i n ag reement with the examples of the Beta and Inverse-Beta 
polymers studied in the text, to analytically continue (Bl) to ( ni,ri 2 ) = (si,S2) € C 2 . Obviously there is an infinite 
number of possible analytical continuations and here we only study arguably the most natural ones (using Euler 
inversion formula r(x)r(l — x) = Tr/sin(Trx ) and (—1)™ = sin(7r(a: + n))/ sin(7rx) for integer n), which we now 
enumerate. 


First type, (ei,e 2 ) = (1,1): 


r(ck + /3) r(a + si)r(/3 + S 2 ) 
r(ck)r(/3) r(a -f /3 -f si -f S 2 ) 


(B2) 


• Second type (£1,62) = (1, — 1): 

_____ = r(l-q) T (1 — a — (3 — si — s 2 )r(/3 + s 2 ) 

U V ~ T(1 - a - p)T(P) F(1 -a- Sl ) 


(B3) 
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• Third type (ei,e 2 ) = (-1, -1): 

= r(a)r(fflr(i-„-<i) r(a + Sl)r( ' 3 + ® )r(1 -<*- H- « -*>) ■ 

• Fourth type, (ei,e 2 ) = (-1,1): 

_ r(a + /3)r(l — a) r (/? + s 2 ) 

r(/3) r(l — a. — si)r(a + p + Si + s 2 ) 


(B4) 


(B5) 


• Fifth type (ei,e 2 ) = (-1,-1): 


u s i v s 2 = r(a + /3)r(i 


a)r(l-/3) 


1 

F(1 - a - si)r(l - /3 - s 2 )r(a + /3 + si + s 2 ) 


• Sixth type (ei, e 2 ) = (1,1): 

r(l - a)r(l -/3) r(l-a-/3- gl - S2 ) 

r(i — a — /3) T(1 - a - si)r(l - P - s 2 ) 


(B 6 ) 


(B7) 


a. First type 


Let us first consider the first type of analytical continuation (B2). Assuming it to be valid on the full complex 
plane, the distribution p(u, v ) can be directly obtained as an ILT 


p(u,v) = 


1 


N, 


a,0 JC i 


ds i 
2l7T 


’c 2 


ds2 U -1-S lv -1-S2 r(q + si)r(/3 + s 2 ) 
2j7t T(a + P + Si + s 2 ) 


(B 8 ) 


where N a ^ a normalization factor, and different contours Ci can be considered. Here we first consider the 

most natural choices: vertical lines passing through the right of all the poles of the integrands located at Si = —a —mi 
and s 2 = —p — m 2 with (mi,m 2 ) € N 2 , e.g. si = —a + 1 + iy\ and s 2 = —/3 + 1 + *y 2 with (yi,y 2 ) € K 2 . We first 
consider the integration on s 2 . If v > 1, the contour can be closed to the right, giving 0 as a result. On the other 
hand, if v < 1, the contour can only be closed to the left and all the poles of s 2 contribute. This shows 


p(u, v) = 6(0 < v < 1 ) 


1 

Nk 




2*7T 


u * y] v 


-1+/3 +m 2 


' at.,3 JC-i - n 

1 m 2=0 

In this expression, one recognizes the taylor expansion 

F(??+ 1 ) 


(-l)™ 2 r(q + s 1 ) 

m 2 ! T(a + si — m 2 ) 


(l-v)* = £(-«) 


fc=0 F (1 + k)T(r] — k + 1) ^ kl T( V -k + 1 ) ’ 


= £ 


(-u) fe r(i? + i) 


with p = a + si — 1, and the convergence is here assured by the fact that v < 1. We thus get 

p(u,v) = 6(0 < v < 1 )-^— [ —u-'-s'v-^n - v )“+ S1 - 1 . 

N a> p Jc 1 2 * 7 T 

And one now recognizes the integral representation of the Dirac 5 distribution: J c ^ w Sl = S(w — 1), and 

p(u,v) = 6(0 < v < 1)—-— u~ 1 v~ 1+ ^( 1 — —— — 1) 


(B9) 


(BIO) 


(BH) 


N, 


a,/3 


= 0(0 < u < + « - 1 ) • 

r (a) r (p) 


(B12) 


Which is exactly the Beta distribution of the Beta polymer, and the normalizibility condition imposes a > 0 and 

P > 0. 
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b. Second type 


Introducing as in the main text 7 = 1 — a — f3, we now study [B3j which reads in these variables: 


u Sl v S2 


We follow the same step as before 


r(7 + (3) r (7 - si - s 2 )T(/3 + s 2 ) 
r (7)r(/3) r(7 + /3-si) 


(B13) 


p{u,v) 


1 f dsi f ds 2 u ^_ aHrl _ S2 r (7 - Si - s 2 )r(/3 + s 2 ) 

Ny,p Jc 1 Jc 2 2i7T r (7 + /? - Si) 


(B14) 


where ^ ■ The contours C\ and C 2 are chosen so as to avoid the poles si = 7 — S 2 +nii and S 2 = — /3 — m 2 

for (mi,m 2 ) £ N 2 , e.g. we choose s 2 = —/? + 1 + iy 2 and Si = y + f3 — 2 + iy 2 with ( 2 / 1 , 2 / 2 ) £ R 2 - Integrating first on 
si and following the same steps as before, we now obtain 


p(u,v) 


= 9(1 < u) 


1 


AL 


, OO 

y 

2in 


7,(3 C 2 “ " mp —0 

0(1 <M) f dS2 „,-l_s 2 „,-l- 7 +s 2 


(~ir 

mi! 


W 


7,/3 J C-2 


2Z7T 


0 (! < u) 1 _i_ 


1 , 


(1 _ a' 


— 1— 7 +S 2 —mi 


/3+s 2 -1 


T(/3 + s 2 ) 

T(/3 + s 2 - mi) 


V u 


Ny,f} 

9(1 < u) 1 


( l --)^5 -(!--)-1 


W 


7,/S 


7(1_ _)^-lj( w _l +u ) 
U 


We thus obtain the same distribution as before, which is indeed normalizable for 7 > 0 and (3 > 0. Note the 
interesting fact that, though the moments of the distribution only exist for ni + n 2 < 7 , the analytical continuation 
of the moments offered by the Gamma function allows us to retrieve the distribution. Though non-rigorous it gives 
some insight to understand why the replica method developed in this paper to retrieve the PDF of the partition sum 
of the associated polymer model works. This is also in agreement with Appendix [D] 


c. Third type 


For the third type (B41, it seems difficult to compute the involved integrals in full generality since they depend on 


the precise position of the poles. However, we now directly exhibit some examples that define proper distributions for 
1 — (a + f3) = 7 > 0. We take u = uw, v = vw with ft, v and w independent random variables distributed with PDF 

1 

1X7)' 


p w (w) = 


Pu(u) = 


Pv(v) = 


1 


r(a) 


T(/3) 


L _ “j k 

Tp) 


,-l+P( p - u 


k\ 

k =0 

L-/ 3 J k 

k =0 


(B15) 


Where [()J denotes the integer part, and the sum appearing in p u (resp. p v ) is present only if a < 0 (resp. /? < 0) 
and regularizes the eventual divergences near the origin. These distributions are singular and only have a few integer 
moments, but their complex moments u Sl v S2 do exist on a domain Re(s\ 4 - S 2 ) < 7 , supplemented by the condition 
|f?e(si)| < 1/2 (resp. |i?e(s 2 )| < 1/2) if a < 0 (resp. /3 < 0), and are there given by (B4). As in the log-Gamma 


and Inverse-Beta cases, these moments can be analytically continued to the full complex plane, opening a way for 
a Bethe ansatz solution of this kind of model. In terms of contours integrals, p u can be obtained using the same 
technique as before with u s = F(o: + si)/r(a), but always choosing a contour of integration as a vertical line passing 
by the origin (and eventually separating the poles of the integrand). It would be of great interests to understand if 
one can obtain exact results for a polymer model defined with these types of weights (e.g. the PDF of log Z t (x)) using 
analytical continuations of other known results. This is left for future work. Notice that these models could well be 
good candidates to fill the void left in the down-left quarter of Fig. [2j 
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d. Other types 


For the other types, one intuitively see that they “lack of poles” in the complex plane (si,S 2 ) £ C 2 to obtain a 


mea ningful result after Laplace inversion, and the corresponding integrals diverge. For the fourth (B5) and fifth cases 
(B 6 ), another argument goes in the same direction. Writing schematically u Sl v S2 = /(si,S 2 ), we have 


7T-yr c d 2 ,, 

( ln «) = ^2/| sl = S2 =o 


7T-yr c 9 2 ,, 

( ln «) = ^|/| S1 = S2 =0 


(B16) 


where (.) 2 denotes the variance. Applying this formula on (B5) and (B 6 ) always leads to negative results which is 


incompatible with In u having a PDF with a second moment. We do not consider here the possibility of such very 
singular distributions. 

For the sixth type (B7), there remains small windows of parameters for which both (Init ) 2 and (Inn ) 2 are positive 


simultaneously, so that this argument is inconclusive. We do not investigate further the possibility of the existence of 
another integrable model here. 


Appendix C: Effect of an additional inverse Gamma weight at the starting point 


Consider the two partition sums, one, noted Z t (x ) and studied in the text, in presence of the additional inverse 
Gamma random variable Wqo on the site x = t = 0, and the other one, Z t (x), in absence of such a weight (which in a 
sense is the true point to point problem). Clearly one has: 

Z t (x) = w 00 Z t {x) (Cl) 


for any x,t where woo and Z t (x) are uncorrelated random variables. 

There are various ways to express one problem into the other. Let us use here the shorthand notation Z = Z t {x) 
and Z = Z t (x). The moments are related as: 




r(7~ s) 

r(7) 




(C2) 


And the Laplace transforms as: 


g—uZ _ £—UWqqZ 


2 

r(7) 


(iiZ)^/ 2 AT 7 (2V / uI) 


(C3) 


Since the l.li.s. is known explicitly as a Fredholm determinant, we see that to obtain P(Z) one needs to invert a 
modified type of Laplace transform involving Bessel functions. 

There is also a useful relation between the CDF’s. Let us define the CDF of In Z, as F{y) = Prob(h\Z < y), and 
the one of In Z. as F(y) = Prob(lnZ < y). Clearly 


F(y) = Prob(\nZ < y) 


rh + dy) 
r(7) 


Hy) 


(Prob(lnZ <y- lniuoo))™™ 


(F(y - lnw 0 o))w 00 = ( e lnw ° ody )w 00 F(y) 


(C4) 

(C5) 


Ltence, knowing F(y) from Fredholm determinants, one can obtain the CFD F(y) as: 

Hy) = , F{y) (C 6 ) 

r (7 + 5 y ) 

an operator which can be interpreted in the sense of a Taylor expansion w.r.t. d y . At fixed 7 in the large time limit 
the rescaling (76 1 renders the term d y smaller by t -1 / 3 . Defining as in (791 


studied in 


IV C 


Fres(z) = Prob(2 3 l0g Z + tCip < z) = F(y = 2 2/3 X lp z - tc v ) 

A(o 


(C7) 


and similarly for F res (z) w.r.t. logZ, we obtain: 

Fres (^) = 


r(7) 


r(7 +2-2/3 A- 1 ^) 


Fres(z) 


(C 8 ) 
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a relation exact for all t, but which for t —► +oo shows that the effect of the operator d y becomes negligible in the 
scaling variable z (we recall that X v ~ t 1 / 3 . The rescaled CDF’s are thus the same, a very intuitive result: the large- 
length is insensitive to such a change in the energy of the first site. The formula allows to calculate the subleading 
corrections. 


In the T = 0 limit the formula (C6) simplifies. Defining Ft= oM = -FVo&((E( tiX ) > r) and similarly for Fp= o(r) in 
presence of the additional exponentially distributed energy random variable at site {x, t) = (0, 0), we obtain from the 
definition (94): 


F T=0 (r) = (1- -d r )F T = 0 (r) 


(C9) 


valid for arbitrary t. At large time the same argument on the rescaled variable z again shows the derivative term to 
be negligible. 

Finally note that such relations have been studied also in the context of stationary models in the KPZ clas d- l^s i 
where it seems also mandatory to add an inverse Gamma variable at the origin in order to obtain a FD representation. 
Its occurence in a point to point problem is, to our knowledge, new. 


Appendix D: Laplace transform Vs Moment Generating function: recall 


In this appendix we briefly recall the idea, discussed iiPS and to which we refer for more details, that leads to the 
conjecture (69). It it best illustrated on the simple problem of obtaining the Laplace transform of ([l]): 


g (X) = e ~ Xu = 


F(7 + /3) 


P+OO 


e -\u _ 


mw) j i 

r(7 + P) f +oc 

mm J i 


due 


— Xu 


,1+7 


du 


^ (—A)" 


n—0 


i - 


,1+7 


1 - 


3-1 


(Dl) 


In this formula, it is obvious that one cannot invert the sum and integrals sign because the different terms converge 
only if n < 7. In which case 


r(7 + ff) 

mw) 


,+°o ( _ A n) 

/ du - :—u 

n! u 


1 


1+7 


1 

1- - 
u 


(—A) n r( 7 + /3)r( 7 — n) 
n\ r (7 + ^-n)r( 7 ) 


(D2) 


Note however that, using the analytical continuation of the Gamma function, the right hand side of (Dl) also makes 
sense for n > 7. We can thus consider the object, called the “moment generating function” defined as 

(-A)” r(7 + /3)r(7-n) 


9n 


>(A) = 


^r(l+n)r (7 + /3-n)r(7) ' 


(D3) 


The question is now to understand how (D3) and (Dl) are related. Let us now rewrite (Dl) using an integral 
representation of the exponential as 


/• + OO 


5(A) = 


du(—l) 


ds 


ds 


X s 


'c 2i sin(7rs) T(1 + s) rt 1+7 


A s 


^*+00 


1 


r(7 + /?) 

mm Ji 

_mmi _. _ 

r(7)r(/3) J e 2 i sin(7rs) T(1 + s) J ± “ u 1+ t 
r ds X s r(7 +/3)F(7 - s) 

Jc 2 i sin(7rs) T(1 + s) r(7 + /3 — s)T(7) 


1 


1 - - 
u 


1 

1- 

u 

3-1 


(D4) 


where here the contour of integration C is a vertical line C = —a + fK with 0 < a < 1. In this way, one can invert 
the different integrals and the results only contains complex moments u s in a region where they are defined. The 
relation between g{ A) and g mom (A) now appears clearly by comparing (D3) and (D4): g( A) can formally be obtained 
by rewriting the sum appearing in g m om{ A) as a Mellin-Barnes transform. Note that when closing the contour of 
integration C on the Re(s) > 0 half-plane in |D4| one obtains two types of poles. A first series coming from the sine 
function that reproduces the series that defines g m om{ A), as well as a second series of terms of the form A 7+rl with 
n £ N coming from the poles of T(y — s ): g( A) is not an analytic function of A. Rewriting the sum appearing in 
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<7mom(A) as a Mellin-Barnes integral thus allows us in some way to retrieve the missing, non-analytic terms that are 
present in the Laplace transform. In the main text we use the same prescription to go from a Fredholm determinant 
formula for g™£ m (u) to a formula for g t , x (u ) by rewriting the sum over m appearing in the expression of the kernel 
(68) as a Mellin-Barnes type integral in (69). Notice that in (68), the sum over m runs from 1 to oo, and the associated 
integral written in (69) is thus chosen as a line that passes through the right of 0 (a trivial modification of the case 
studied here), and to the left of 7 to avoid crossing a pole. 
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